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Abstract. Between June 2006 and September 2009, an in- 
strumented light aircraft measured over 400 vertical profiles 
of aerosol and trace gas properties over eastern and cen- 
tral Illinois. The primary objectives of this program were 
to (1) measure the in situ aerosol properties and determine 
their vertical and temporal variability and (2) relate these air- 
craft measurements to concurrent surface and satellite mea- 
surements. The primary profile location was within 15 km 
of the NOAA/ESRL surface aerosol monitoring station near 
Bondville, Illinois. Identical instruments at the surface and 
on the aircraft ensured that the data from both platforms 
would be directly comparable and permitted a determina- 
tion of how representative surface aerosol properties were 
of the lower column. Aircraft profiles were also conducted 
occasionally at two other nearby locations to increase the 
frequency of A-Train satellite underflights for the purpose 
of comparing in situ and satellite-retrieved aerosol data. 
Measurements of aerosol properties conducted at low rel- 
ative humidity over the Bondville site compare well with 
the analogous surface aerosol data and do not indicate any 
major sampling issues or that the aerosol is radically dif- 
ferent at the surface compared with the lowest flyby alti- 
tude of ~ 240 m above ground level. Statistical analyses of 
the in situ vertical profile data indicate that aerosol light 
scattering and absorption (related to aerosol amount) de- 
creases substantially with increasing altitude. Parameters re- 
lated to the nature of the aerosol (e.g., single-scattering 
albedo. Angstrom exponent, etc.), however, are relatively 
constant throughout the mixed layer, and do not vary as much 
as the aerosol amount throughout the profile. While individ- 


ual profiles often showed more variability, the median in situ 
single-scattering albedo was 0.93-0.95 for all sampled al- 
titudes. Several parameters (e.g., submicrometer scattering 
fraction, hemispheric backscattering fraction, and scattering 
Angstrom exponent) suggest that the fraction of smaller par- 
ticles in the aerosol is larger near the surface than at high 
altitudes. The observed dependence of scattering on size, 
wavelength, angular integration range, and relative humidity, 
together with the spectral dependence of absorption, show 
that the aerosol at higher altitudes is larger, less hygroscopic, 
and more strongly absorbing at shorter wavelengths, suggest- 
ing an increased contribution from dust or organic aerosols. 
The aerosol profiles show significant differences among sea- 
sons. The largest amounts of aerosol (as determined by me- 
dian light extinction profile measurements) throughout most 
of the sampled column were observed during summer, with 
the lowest amounts in the winter and intermediate values 
in the spring and fall. The highest three profile levels (3.1, 
3.7, 4.6 km), however, showed larger median extinction val- 
ues in the spring, which could reflect long-range transport 
of dust or smoke aerosols. The aerosols in the mixed layer 
were darkest (i.e., lowest single-scattering albedo) in the fall, 
in agreement with surface measurements at Bondville and 
other continental sites in the US. In situ profiles of aerosol 
radiative forcing efficiency showed little seasonal or verti- 
cal variability. Underflights of the CALIPSO satellite show 
reasonable agreement in a majority of retrieved profiles be- 
tween aircraft-measured extinction at 532 nm (adjusted to 
ambient relative humidity) and CALIPSO-retrieved extinc- 
tion, and suggest that routine aircraft profiling programs can 
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be used to better understand and validate satellite retrieval 
algorithms. CALIPSO tended to overestimate the aerosol ex- 
tinction at this location in some boundary layer flight seg- 
ments when scattered or broken clouds were present, which 
could be related to problems with CALIPSO cloud screening 
methods. The in situ aircraft-collected aerosol data suggest 
extinction thresholds for the likelihood of aerosol layers be- 
ing detected by the CALIOP lidar. In this study, aerosol lay- 
ers with light extinction (532 nm) values >50 Mm -1 were 
detected by CALIPSO ~ 95 % of the time, while aerosol lay- 
ers with extinction values lower than 10 Mm -1 had a detec- 
tion efficiency of <2%. For all collocated comparison cases, 
a 50 % probability of detection falls at an in situ extinction 
level of 20-25 Mm -1 . These statistical data offer guidance 
as to the likelihood of CALIPSO’s ability to retrieve aerosol 
extinction at various locations around the globe. 


1 Introduction 

Unlike the major long-lived trace gases, aerosols are not dis- 
tributed uniformly in the atmosphere. Locations downwind 
of major particle sources such as deserts, biomass burning re- 
gions and large cities often have heavy atmospheric aerosol 
burdens while areas remote from major sources have rela- 
tively low aerosol loadings. The composition and size dis- 
tribution of these atmospheric particles, and thus their opti- 
cal and hygroscopic properties, also vary widely from place 
to place and over time. These aerosol properties depend not 
only on the source emissions but also the subsequent atmo- 
spheric processing including condensation, coagulation and 
removal, which can be quite different for different aerosol 
types. This inhomogeneity in aerosol amount and charac- 
ter coupled with the relative paucity of sustained observa- 
tions around the globe makes it difficult to estimate the di- 
rect aerosol radiative forcing effect on global climate (e.g., 
Charlson et al., 1991; Yu et al., 2006). 

In order to better understand aerosol radiative forcing 
and the effects of aerosols on the global and regional cli- 
mate, long-term measurements of aerosol optical properties 
are being made by the National Oceanic and Atmospheric 
Administration (NOAA) Earth System Research Laboratory 
(ESRL) (e.g., Bodhaine, 1983, 1995; Bodhaine and Dutton, 
1993; Delene and Ogren, 2002) and other organizations (e.g., 
the Global Atmosphere Watch network coordinated by the 
World Meteorological Organization) at locations around the 
world. Many major surface regions remain undersampled, 
however, and very few long-term measurement efforts have 
been made above the surface. Despite these dedicated aerosol 
monitoring programs providing important model initializa- 
tion and validation data, the largest uncertainties in mod- 
eling climate change remain those associated with aerosols 
(e.g., IPCC, 2007; Myhre, 2009). Part of the problem is that 
it is unknown to what extent the surface measurements at 


many locations represent the aerosols above them in the tro- 
posphere. For reliable aerosol radiative forcing estimates to 
be made, knowledge of the aerosol optical properties, at least 
up to the top of the mixed layer where most of the aerosols 
reside, is necessary. 

Remote sensing methods, including aerosol optical depth 
retrievals and lidar measurements from both satellite- and 
surface-based platforms, are clearly useful in helping to un- 
derstand the vertical distributions of aerosols. These methods 
have their limitations however. While polar-orbiting satellites 
provide broad spatial coverage, their sensors require careful 
validation before meaningful information on the nature of the 
aerosols can be extracted from the retrievals. It has also re- 
cently been shown that errors in the assumed aerosol pro- 
files can cause significant errors in the aerosol optical thick- 
ness retrieved from satellites (Rozwadowska, 2007). Inver- 
sion algorithms utilizing remote sensing measurements from 
ground-based sun-sky radiometers (e.g., multi-wavelength 
Cimel sun photometers) are promising (Holben et al., 1998, 
2001; Eck et al., 2003, 2005), but spatial coverage is quite 
limited compared to satellite measurements. Additionally, 
the aerosol properties obtained from sun photometer mea- 
surements are representative of the entire column rather than 
those at any specific altitude, have difficulties with the re- 
trievals of intensive properties for lower aerosol loadings, 
and are normally limited to measure in clear sky conditions 
(e.g., Dubovik and King, 2000). 

Lidar routinely provides profiles of aerosol extinction and 
backscattering (“extensive” aerosol properties, related to the 
amount of aerosol present, Ogren, 1995). At this time, how- 
ever, lidar retrievals of the aerosol intensive properties (re- 
lated to the nature of the aerosols rather than the amount, 
Ogren, 1995) necessary to determine radiative forcing re- 
main less certain (Muller et al., 2001; Veselovskii et al., 
2005). A large data set has in the last few years been col- 
lected by NASA using their High Spectral Resolution Li- 
dar (HSRL) instrument and has shown promise in measuring 
some intensive aerosol properties including lidar ratio, depo- 
larization, and backscatter wavelength dependence (Rogers 
et al., 2009). 

Instrumented aircraft are capable of making many of the 
in situ aerosol measurements necessary for validating the 
remote sensing methods. In addition, they provide valuable 
data for initialization and testing of global aerosol transport 
models and for relating surface measurements to those in 
the overlying column. Most of the historical airborne mea- 
surements have been during short “deployment-based" stud- 
ies (e.g., Verver et al., 2000; Ramanathan et al., 2001; Fer- 
rare et al., 2006; Shinozuka et al., 2007) intensely study- 
ing the atmosphere for relatively brief periods of time (e.g., 
several weeks). The duration of these projects typically lim- 
its their usefulness in comparing the airborne data with sur- 
face or remote sensing climatologies and is insufficient to 
determine seasonal aerosol variability or long-term trends 
to assess climatic effects. Unfortunately, not many aircraft 
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programs have been conducted even intermittently in one re- 
gion over multiple years (e.g., Taubman et al., 2006; Hains et 
al., 2008), and even fewer can be considered true long-term 
monitoring programs (e.g., Andrews et al., 2004, 2011). 

In an effort to characterize when, how often, and under 
what conditions surface aerosol measurements are represen- 
tative of the lower atmospheric column, ESRL recently con- 
ducted long-term aircraft measurement programs over two 
heavily-instrumented surface aerosol monitoring stations. 
These sites were the Department of Energy /Atmospheric Ra- 
diation Measurement (DOE/ARM) program Southern Great 
Plains (SGP) Central Facility near Lamont, Oklahoma, and 
the NOAA surface aerosol monitoring station (BND) near 
Bondville, Illinois. The primary objectives of these programs 
were to obtain statistically robust data sets of the vertical dis- 
tribution of aerosol properties for use by modelers in the 
evaluation of aerosol radiative forcing and to relate these 
properties to those measured by similar or identical instru- 
ments at the surface. To accomplish these primary goals, ded- 
icated aerosol systems were built into small Cessna aircraft 
and were designed using identical inlets and many of the 
same instruments, so that measurements would be directly 
comparable between the two aircraft systems and between 
surface and airborne stations. The Oklahoma aircraft flew 
597 complete research profiles between March 2000 and De- 
cember 2007, and results have been reported in Andrews et 
al. (2004) and Andrews et al. (2011). The Illinois aircraft is 
the NOAA/ESRL Airborne Aerosol Observatory (AAO), and 
it flew 2-3 times per week between lune 2006 and September 
2009 for a total of 401 research profiles. 

The measurements taken aboard the two aircraft were used 
to answer the following scientific questions: 

- How do aerosol properties in the atmosphere over a ru- 
ral region vary over different time scales (e.g., monthly 
or seasonally)? 

- How do aerosol properties change with altitude in the 
lower atmospheric column? 

- At what times, and under what conditions can surface- 
based measurements of these properties be used to cal- 
culate the direct aerosol radiative forcing from a mea- 
sured aerosol optical depth? 

A secondary objective of the AAO project was to contribute 
to the verification of aerosol remote sensing retrieval algo- 
rithms used by several satellites, including the “A-Train” 
AQUA and CALIPSO satellites and the TERRA satellite. 
This was accomplished by coordination of aircraft activi- 
ties with satellite overflights and generation of value-added 
products for comparison with the lidar and column-average 
sunphotometer data. The satellite measurements require “air 
truth” in situ measurements to identify and constrain sys- 
tematic errors in the retrievals of aerosol physical properties. 
Improvements in these algorithms will improve the quality 


of height dependent information from CALIPSO and other 
backscatter lidars. 

This paper presents the statistics of the > 3-yr AAO record 
of in situ aerosol optical property measurements. The verti- 
cal profiles of aerosol properties over the BND site are useful 
in determining the statistical variability of aerosol properties 
with time and altitude in central Illinois. While the AAO data 
are not of sufficient duration to determine long-term aerosol 
trends, they are useful in assessing seasonal variability and 
also in relating the atmospheric and surface aerosol mea- 
surements. Measurements of in situ AAO light extinction are 
compared with CALIPSO satellite lidar (CALIOP) extinc- 
tion data, and the results suggest a cost-effective way to ver- 
ify satellite-borne lidar retrievals. Comparisons of AAO light 
extinction with aerosol optical thickness measurements from 
the Bondville AERONET sun photometer are reported in an- 
other paper (Esteve et al., 2012). 


2 Methods 

The in situ aerosol data were obtained by flying an in- 
strumented light aircraft (Cessna T206H) near NOAA’s re- 
gional aerosol monitoring station near Bondville, Illinois, 
(40.053° N, 88.372° W, 220m above mean sea level - a.s.l.). 
This site is an anthropogenically perturbed, continental sta- 
tion located at the Illinois State Water Survey’s Bondville En- 
vironmental and Atmospheric Research Site (BEARS). It is 
located 6.5 km south of Bondville (population ~ 450), 16 km 
southwest of Champaign-Urbana (population ~ 230 000), 
and is surrounded by corn and soybean fields. The prevail- 
ing wind directions over an entire year at the site cover 
a range from S to WNW (i.e., predominantly upwind of 
Champaign-Urbana). Climatologies of surface aerosol opti- 
cal properties observed at the Bondville site have been re- 
ported by Koloutsou-Vakakis et al. (2001) and Delene and 
Ogren (2002). 

A schematic of instrument rack locations in the Cessna 
aircraft is shown in Fig. 1. The locations of instruments are 
reported in Table 1, which also lists details about the mea- 
surements and the species or parameter measured. The main 
aerosol sample line running to the optics rack (O) and the 
ABC rack system is shown as a red line. 

The aerosol sample inlet was made of anodized aircraft- 
grade Al alloy, based on the shrouded inlet design of the 
University of Hawaii (Clarke et al., 2004; McNaughton et al., 
2007). The inlet was positioned on the starboard wing about 
2 m from the fuselage, well outside of the propeller wash. 
The exit plane of the inlet was modified for a smooth tran- 
sition to 22.2 mm internal diameter (ID) stainless steel tub- 
ing, and this inlet line entered the leading edge of the wing. 
The front of the aerosol inlet extended ~ 50 cm forward of 
the leading edge of the wing, sufficient to get into the free air 
stream at ~ 50 m s -1 airspeed with little disturbance from the 
forward projection of the wing wake based on tufting studies 
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Table 1 . Major instruments and primary measurements on the AAO aircraft. 


Aircraft Location 

Instruments 

Measurements/Species 

Wavelengths (nm) 

Rack O 

3-A. Nephelometer 

cr S p and cr^sp ( no s * ze cu f l° w RH) 

450, 550, 700 


Three 1-A Nephelometers 

OsplUp < 1 pm, 3 different RHs) 

545 

Rack A 

CPC 

Total particle number (Dp > lOnm) 

- 


SEMS 

Particle size distribution (22 nm < Dp < 480 nm) 

- 

Rack B 

PILS 

Soluble fraction major ions 

- 

Rack C 

PSAP 

(r a p (no size cut) 

467, 530, 660 


Ozone monitor 

Ozone mixing ratio 

- 

Rack F 

Gas flask sampler 

Trace gases 

- 



Fig. 1. Schematic of the AAO Cessna T206H aircraft showing 
aerosol inlet and instrument racks. Main aerosol sample line is 
shown in red. Lettered areas are described in the text and in Table 1 . 


done at that location. Those same tufting studies guided our 
angling of the inlet downward by several degrees relative to 
the angle of the wing to match the orientation of the stream- 
lines in level flight when the aircraft was burdened with a 
low-to-medium load of fuel. Adjustment of the inlet orien- 
tation during flight was not possible, so non-axial sampling 
by up to a few degrees is possible early in the flight. Sam- 
pling later in the profiles, especially that conducted down in 
the boundary layer where most of the aerosol resides, was 
expected to be close to isoaxial. 

The aerosol sample flow rate was maintained at ~ 60 1pm 
so that the inlet would sample isokinetically at the nominal 
50 ms -1 research speed of the Cessna. Upon entering the 
wing the inlet tubing made a 90° gentle bend (7.6-cm radius 
of curvature) toward the fuselage followed by an approxi- 
mately level run to the cabin. Inside the cabin, the sample 
flow encountered a flow splitter that directed a portion of the 


flow (16.7 1pm) downward into the optics rack. The flow rate 
was controlled by a volume flow controller (Brechtel Man- 
ufacturing Inc.), and this flow control was necessary due to 
a downstream aerosol cyclone. The remainder of the aerosol 
flow was passed across the cabin ceiling to the ABC rack 
system. 

The performance of the University of Hawaii (UH) in- 
let has been described in detail (Huebert et al., 2004; Mc- 
Naughton et al., 2007). The results indicate a 50 % passing 
efficiency for ~ 5 pm aerodynamic diameter particles for the 
UH inlet, including its short length of inlet tubing, at an air- 
speed of 120 ms~’. Scaling this result to our lower airspeed 
suggests a 50 % passing efficiency through this shrouded dif- 
fuser inlet for particles of 6-7 pm aerodynamic diameter. The 
particle passing efficiency of the AAO sample line tubing 
downstream of the inlet was not determined experimentally, 
but was calculated using the widely available aerosol calcula- 
tor spreadsheet from Paul Baron (e.g., http://aerosols.wustl. 
edu/AAARworkshop08/html/calculator.htm). For this exer- 
cise, we neglected electrostatic (all tubing was conductive) 
and thermophoretic losses, and calculated the combined ef- 
fect of gravitational, inertial, and diffusional losses in the tur- 
bulent flow conditions inside the AAO inlet tubing for parti- 
cles of different aerodynamic diameters. Figure 2 indicates 
that the inlet plumbing to the inlet of the reference neph- 
elometer in the optics rack efficiently passed particles over 
the range of ~ 20 nm up to ~ 5 pm. The particle transmission 
efficiency to the flow splitter in the ABC rack system experi- 
enced only slightly higher transmission losses at the low and 
high ends of the particle size distribution due to the longer 
run of tubing. The combination of literature results and loss 
calculations indicates that the AAO inlet system (including 
the diffuser inlet plus sample line tubing) efficiently passes 
particles up to about 5 pm aerodynamic diameter. 

The optics rack was positioned where the co-pilot’s seat 
would normally be and housed a three -wavelength (3-A.) in- 
tegrating nephelometer (TSI Model 3563) for measurement 
of the total light scattering coefficients (er sp ) and hemispheric 
backscattering coefficients (ffbsp)- Measurement details and 
uncertainties for the TSI nephelometer have been described 
elsewhere (Anderson et al., 1999; Sheridan et al., 2002). 
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Fig. 2. Particle transmission efficiency through the aerosol sample 
lines between the main inlet and the Optics and ABC racks (cal- 
culated using Paul A. Baron’s aerosol calculator spreadsheet, AE- 
ROCALC.XLS). Curves represent the combined effects of gravi- 
tational, inertial and diffusional losses in a turbulent flow environ- 
ment. 

Ambient aerosol properties for this program (e.g., single- 
scattering albedo ( ojq ) > 0.7 and minor coarse mode frac- 
tions, as discussed by Massoli et al., 2009) were such that 
the routine nephelometer correction (i.e., truncation) meth- 
ods described in Anderson and Ogren (1998) were appropri- 
ate and, thus, were applied in this study. 

Sample aerosols entering the optics rack passed through 
the 3-/, nephelometer (controlled at a relative humidity (RH) 

< 40 % by a small heater) and then entered a sharp cut cy- 
clone (BGI Inc. model SCC 2.229-PM1), which produced 
a 1 -pm aerodynamic diameter size cut aerosol at a volume 
flow rate of 16.71pm. This aerosol was then fed in parallel 
into three single- wavelength (1-7.) integrating nephelometers 
(Radiance Research Model M903), each operating at a vol- 
ume flow rate of ~ 5.61pm, for measurement of the submi- 
crometer CT sp . The 1-3. nephelometer data were also corrected 
for truncation using the correction factors presented in An- 
derson and Ogren (1998). The question of whether the An- 
derson and Ogren (1998) corrections, which were published 
to correct TSI nephelometer truncation errors, can also be 
used for the Radiance Research nephelometers has been an- 
swered in the recent study of Muller et al. (2009). This study 
shows that the measured illumination functions for the Radi- 
ance Research and TSI nephelometers are very similar, and 
these coupled with the internal instrument geometries pro- 
vide truncation correction factors for the Radiance Research 
nephelometer that agree quite closely (within a few percent 
up to 1-pm aerodynamic diameter) with the original Ander- 
son and Ogren (1998) corrections. The three 1-7. nephelome- 
ters were each held at a different stable RH through the use of 
small proportional-integral-derivative (PID) controllers (low: 

< 40 % RH, med: 65 % RH, high: 85 % RH) for determina- 
tion of the submicrometer aerosol hygroscopic growth factor. 


The remainder of the 601pm inlet flow was passed across 
the cabin to the ABC rack system. Rack A housed a conden- 
sation particle counter (CPC, TSI Model 3010) for measur- 
ing condensation nucleus (CN) concentration and a scanning 
electrical mobility sizer (Brechtel Manufacturing Inc. Model 
SEMS) for measuring aerosol size distributions for parti- 
cles between 22 and 480 nm. Rack B contained a Particle - 
Into-Liquid Sampler (PILS, Brechtel Manufacturing Inc.) 
for post-flight determination of aerosol ionic composition. 
Aerosol samples of ~ 3 min duration were collected by in- 
jection into small sealed vials at 5 min intervals and sub- 
sequently sent to NOAA’s Pacific Marine Environmental 
Laboratory for major ion analysis. Rack C housed a 3-7. 
filter-based light absorption instrument (Radiance Research 
Model PSAP) for measurement of the aerosol light absorp- 
tion coefficient (er ap ) and a continuous ozone monitor (2B 
Technologies Model 205). Rack F was a removable suitcase 
sampler (http://www.esrl.noaa.gov/gmd/ccgg/aircraft/index. 
html) that contained 12 glass flasks for trace gas collec- 
tion and post-flight analysis at NOAA/ESRL (e.g., CO 2 , CO, 
CH 4 , N 2 O, H 2 , SF(,, isotopes of CO 2 and CH 4 , multiple halo- 
and hydrocarbons). A gas sample was collected near the mid- 
point of each of the 10 level segments of each research pro- 
file. Typically, a second flask would be collected on one of 
the higher altitude runs to assess analytical precision and also 
on the lowest leg near the BND site on the way back to the 
airport. The gas inlets for the ozone monitor and the suitcase 
sampler were rear-facing inlets located on a modified win- 
dow plate below the port wing. 

PSAP measurements were corrected for sample area, flow 
rate and nonidealities in the manufacturer’s calibration (Bond 
et al., 1999; Ogren, 2010). Measurement uncertainties for 
the PSAP have been described in detail elsewhere (Anderson 
et al., 1999; Bond et al., 1999; Sheridan et al., 2002). The 
AAO PSAP was modified to accommodate a small (~ 5W) 
heater on its internal inlet line just upstream of the optical 
block. The heating was gentle, typically only a few degrees, 
but the metal block was kept warm relative to the sample so 
the RH remained low. Our own tests suggest that the heater 
kept the RH at the filter below 40 % most of the time and an 
RH of 50 % at the filter would have been exceeded only dur- 
ing sampling of very humid air, which was not common in 
this program. The heater was not actively controlled to main- 
tain a specific low RH, but RH variability at low RH is not 
believed to influence the measurements as strongly as RH 
variability at high RH. Figure 2 shows the calculated par- 
ticle transmission efficiency through the entire sample line 
to the PSAP inlet. This calculation shows higher particle 
losses at the ends of the size distribution for the PSAP rel- 
ative to the ABC rack system due to an added length of 
smaller diameter (~ 6 mm ID) conductive tubing between the 
ABC splitter and the PSAP inlet. The PSAP inlet line parti- 
cle passing efficiency compares favorably, however, to the 
transmission efficiency of the inlet tubing to the optics rack 
over the optically-important 20 nm to 2 pm size range. At 


www.atmos-chem-phys.net/12/11695/2012/ 


Atmos. Chem. Phys., 12, 11695-11721, 2012 


11700 


P. J. Sheridan et al.: Vertical profiles of aerosol optical properties over central Illinois 


aerodynamic diameters of 20 nm and 2 pm, > 95 % of all 
sampled particles reached the AAO instmments. This is im- 
portant for the determination of the single-scattering albedo, 
<uo, which has both a scattering and an absorption compo- 
nent. Inlet losses should be similar to those experienced by 
the nephelometer for valid comparisons between scattering 
and absorption. 

Unfortunately the SEMS and PILS measurements were 
quite intermittent in nature due to numerous problems as- 
sociated with operating in an aircraft environment. Vibration 
appears to have been one of the major issues, causing inter- 
nal circuit boards and/or connectors in the SEMS to come 
loose during flights and the injection needle for the PILS to 
frequently stick in the retracted position. Because we have 
such a limited subset of the SEMS and PILS data available, 
only aerosol optical properties and particle number concen- 
trations are included in the long term statistics presented in 
this paper. 

The basic AAO aerosol optical measurements were used 
to derive a number of other key aerosol properties that are 
used in radiative transfer models and for the determination of 
direct aerosol radiative forcing (Haywood and Shine, 1995). 
These calculated properties and equations are listed in Ta- 
ble 2. The aerosol light extinction coefficient, cr e p, is the 
sum of the scattering and absorption coefficients at any given 
wavelength. The scattering and absorption Angstrom expo- 
nents (Ensor et al., 1972; Bodhaine, 1983; Moosmiiller et 
al, 2009), a s and a a , respectively, describe the wavelength- 
dependence of light scattering and absorption based on a 
power law relationship. For comparison with historical data, 
we used a common two-point power-law expression for cal- 
culating the a s . We used a more robust 3-point fit (Virkkula 
et al., 2011) for the a a determination. The submicrometer 
scattering ratio, 7? sp , is the fraction of light scattering at a 
given wavelength due to aerosol particles smaller than 1 pm 
aerodynamic diameter. As discussed above the inlet and in- 
let tubing efficiently pass particles up to ~ 5 pm aerodynamic 
diameter, so the “total” size fraction for AAO measurements 
includes particles up to ~ 5 pm. The a>o is the fraction of light 
extinction due to scattering. 

A number of other parameters describe the angular scat- 
tering properties of the aerosol. The <7b sp and (T sp measure- 
ments are used to calculate the fraction of light scattered in 
the backward hemisphere (the hemispheric backscatter frac- 
tion, b) at each wavelength. The hemispheric backscatter 
fraction is related to the aerosol asymmetry parameter, g, 
which describes the angular scattering phase function, and 
is estimated from the empirical relationship presented in An- 
drews et al. (2006). The upscattered fraction /, or the frac- 
tion of incoming light scattered backward to space by at- 
mospheric aerosols, has been related to b by Wiscombe and 
Grams (1976), and a second-order curve fit of the points in 
their Fig. 3 as reported in Sheridan and Ogren (1999) pro- 
vides the parameterization shown in Table 2. 



Fig. 3. Map of the study area showing the three AAO profile lo- 
cations, the approximate A-Train satellite tracks, the BND surface 
monitoring site, and the AAO base of operations (CMI airport). 


The submicrometer aerosol scattering hygroscopic growth 
factor, /(RH), for each altitude was estimated using a subset 
of measurements from the three 1 — X (545 nm) nephelome- 
ters. These nephelometers measured submicrometer scatter- 
ing coefficients at three fixed RH’s (< 40 %, ~ 65 %, and 
~ 85 %), and the mean scattering growth function for each 
flight segment was determined by a power-law fit of the 
three sets of measurements. The /(RH) equation in Table 2 
shows how the fit parameter y is used to derive /(RH), 
which by definition in this work is the relationship between 
°s P ,rh=85% relative to cr sp ,R H = 40 %- In this equation RH wet = 
85 and RHd ry = 40. The y parameter describes the curva- 
ture of the RH-dependence of scattering (Kasten, 1969) and 
has been used to characterize aerosol hygroscopicity in other 
studies (e.g., Gasso et al., 2000; Sheridan et al., 2002; An- 
drews et al., 2004; Hains et al., 2008). The same equation 
is used to adjust the submicrometer aerosol scattering values 
measured on each flight segment in the low-RH, 1-k neph- 
elometer to ambient RH conditions. For this adjustment, the 
median y value for each flight level acquired for the sub- 
set of flight segments with valid hygroscopic growth mea- 
surements (i.e., all three nephelometers functioning properly) 
was used with R H j rv now being the RH measured in the 
low-RH nephelometer, RH we t being the ambient RH, crsp,40% 
being replaced by the er sp d ry (he., °sp measured in the low- 
RH nephelometer) and solving for the cr sPiWe t (ambient er sp ). 
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Table 2. Parameters and equations used for calculation of aerosol optical properties. 


Parameter 

Equation 

Extinction coefficient 

"ep — "sp + "ap 

Scattering Angstrom exponent 

a s = — log(cr S px t / cr s P ,A 2 )/ l°g(7 1 / k2) 

Absorption Angstrom exponent 

a a = -log(cr apA ( ( )/log(k(D) 

Submicron scattering ratio 

R sp — "sp.D p <l pm/"sp, total 

Hemispheric backscatter fraction 

h = "bsp / "sp 

Asymmetry parameter 

g = - 7 . 1439 b 3 + 7.4644Z> 2 - 3.9636b + 0.9893 

Upscattered fraction 

/3 = 0.08 17 + 1 .84956 - 2.9682 b 2 

Single-scattering albedo 

o> 0 — "sp/ "ep = "sp/ ("sp + "ap) 

Hygroscopic growth factor 

/( RH) = "sp,85 %/"sp,40% = ((1 - (RHwet/100))/(l - (RH dry /100)))-^ 


Table 3. AAO profile locations and details. 


Profile ID 

Profile Type 

Location (Lat., Lon.) 

Number of Profiles 

LODG 

Routine profile at Lodge location 

40.117° N, 88.567° W 

297 

TERA 

TERRA satellite underflight at Lodge location 

40.117° N, 88.567° W 

41 

AQB1 

A-Train satellite underflight on B Track 

40.717° N, 90.408° W 

20 

AQC1 

A-Train satellite underflight at Lodge location (C Track) 

40.117° N, 88.567° W 

27 

AQD1 

A-Train satellite underflight on D Track 

40.230° N, 87.128° W 

16 


The effect of using the median altitude-specific y value as 
opposed to a y value determined for each flight segment is 
discussed in Sect. 3.2. 

The basic sampling strategy employed for the AAO project 
was to use a dedicated light aircraft to sample atmospheric 
aerosols at predetermined altitudes above an instrumented 
ground station that makes similar or identical measurements 
and/or below a satellite making aerosol measurements. Ta- 
ble 3 provides the locations of the AAO profiles and the num- 
ber of profiles flown at each location, while Fig. 3 shows a 
map of the area with profile locations and A-Train satellite 
tracks indicated. 

To be able to make valid comparisons with the remote col- 
umn sensors, the AAO aircraft had to reach altitudes above 
the vast majority of the aerosols, and unbroken clouds could 
not be present. We have no long-term lidar data from this re- 
gion to assess aerosol extinction vs. altitude; however, sev- 
eral studies suggest that large extinction at high altitudes 
is rare over the US. Yu et al. (2010) reports a summary of 
CALIPSO data over the eastern US and shows that neg- 
ligible extinction was observed at altitudes of 4.6-10 km 
in 2007. Likewise, Vernier et al. (2011) demonstrate that 
the zonally-averaged aerosol optical depth (AOD) of strato- 
spheric aerosols at 20-30 km altitude at 20-50° N latitude 
during the study period was < 0.005. Finally, a detailed 
Raman lidar study was done at another instrumented mid- 
latitude continental site (the DOE/ARM SGP site. Tinner 
et al., 2001) where the AAO sister aircraft operated for 
many years. These lidar data represent a wide variety of 


atmospheric conditions and analyses were analyzed based 
on time of day, season, and integrated amount. The re- 
sults indicate that for AOD < 0.4, our highest flight altitude 
(~4.6kma.s.l.) would be above > 95 % of column aerosol 
extinction. Based on analysis of our AAO profile data, sig- 
nificant aerosol loadings at the highest levels of our profiles 
were rarely observed, so the SGP lidar study results may be 
representative of other mid-continental sites like BND. Inter- 
estingly, the lidar observed some elevated aerosol layers dur- 
ing high AOD (0.6-1 .5) cases; these layers were at and above 
the highest flight altitudes and occurred most frequently in 
the springtime at the SGP site. This observation is consis- 
tent with our springtime in situ aerosol profiles over central 
Illinois, which are reported below. 

Our strategy for operating when clouds were present was 
for the AAO pilot to assess the current extent of cloudiness 
in the region and the cloud forecast prior to takeoff. For rou- 
tine profiles conducted near the small town of Lodge (with 
no satellite overflight), one of two conditions needed to be 
met; either cloud coverage was broken (clouds covering 6/10 
to 9/10 of the sky) or less with a forecast for no increase 
in cloudiness, or there was a shallow solid layer of (non- 
precipitating) cloud that spanned no more than three adja- 
cent flight levels (in which case those levels were not sam- 
pled). For this reason, some AAO profiles are missing 1-3 
level segments. For satellite overflight profiles, cloud frac- 
tion must have been no more than scattered (1/10 to 5/10 sky 
coverage) without a deteriorating forecast. Unfortunately, the 
forecast was not always correct and there were flights with 
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Fig. 4. A typical aircraft vertical profile over the Lodge profile site. Altitudes are relative to mean sea level. A 5-min low pass was typically 
conducted over the Bondville surface monitoring station on the way back to the airport. 


more cloudiness than anticipated, which negatively affected 
the satellite data retrievals and their comparisons with the 
aircraft data. In general, we chose to fly on sunny to partly 
cloudy days; thus our results may be biased toward lower 
ambient RH conditions. 

A diagram depicting a typical research profile over the 
Lodge sampling location is shown in Fig. 4. For all AAO 
profiles, level flight segments centered on the profile loca- 
tion were conducted at 4580, 3660, 3050, 2440, 1830, 1520, 
1220, 920, 610 and 460ma.s.l. (corresponding to 15000, 
12 000, 10000, 8000, 6000, 5000, 4000, 3000, 2000 and 
1500feeta.s.l.). The altitudes flown were pressure altitudes 
based on the aircraft altimeter, which in the US are denoted 
in feet rather than meters. The upper five segments (down to 
and including the 1830 m level) were 10 min in duration to 
give better averaging statistics and reduced uncertainties in 
the lower aerosol-loading portion of the profiles. The lower 
five segments were 5 min in duration. The Lodge site was 
chosen as our primary profile location, rather than directly 
over the BND station, due to FAA concerns over the prox- 
imity to the flight approach path for the Champaign-Urbana 
Willard (CMI) airport. The profile location near Lodge is 
about 15 km west-northwest of BND. On these Lodge pro- 
files (labeled LODG, AQC1, and TERA in Table 3), a low- 
altitude (~460ma.s.l, ~240m above ground level - a.g.l.) 
flyby of the BND site was conducted on the way back to 
the airport to relate the aerosol measurements in the vertical 
profile to those at the surface site. Therefore, all profiles con- 


ducted at the Lodge site have two low-altitude segments; one 
at the bottom of the Lodge profile and one centered ~ 15 km 
away directly over the BND site. 

The aerosol system was designed for semi-automated op- 
eration such that the pilot had only to turn it on for a pre-flight 
warm-up, install gas flasks and expendable supplies such as 
filters and fluids, go through a system checklist to ensure 
that instruments were working properly, and then proceed to 
takeoff. Once airborne, the only interaction of the pilot with 
the research equipment was to press a button near the mid- 
dle of each level-flight segment to initiate collection of a gas 
flask sample. On each flight a nephelometer filtered-air check 
to determine nephelometer background values was automat- 
ically performed just after takeoff and on the way back to 
the airport after completing the profile. Sampling times were 
limited to daylight hours, but otherwise were quite variable, 
so that cumulative data represent all daylight hours, days of 
the week, and seasons of the year. During satellite under- 
flights, takeoff was scheduled such that the airplane was sam- 
pling in the boundary layer when the satellite passed over. 

3 Results 

Figure 5 is a histogram showing the number of each type 
of vertical profile per quarter-year. The AAO program av- 
eraged about 30 research profiles per quarter or about one 
profile every 3 days. An exception was the first quarter of 
2008, where an engine problem kept the aircraft grounded 
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AAO Flight Frequency, June 2006 - September 2009) 


Total Number of Research Flights (excluding Test Flights) = 401 



Fig. 5. AAO flight frequency distribution by quarter-year. The rou- 
tine Lodge and TERRA satellite profiles were all conducted over the 
Lodge profile location. Each A-Train satellite profile was conducted 
at one of the three AAO profile locations listed in Table 3. 


for nearly 2 months. Profiles conducted during overpasses of 
the A-Train constellation and TERRA satellites are shown 
in yellow and red, respectively. About 26 % of the research 
profiles occurred during satellite overpasses, and this frac- 
tion of satellite flights was determined by weather and the 
more stringent cloud fraction requirement, along with pilot 
availability and the satellite overpass schedule. 

One of the objectives of the AAO program was to deter- 
mine statistically how representative the surface measure- 
ments at the BND site were of those in the lower column. In 
Fig. 6, the low-altitude flight segments conducted during fly- 
bys of the BND site at ~ 240 m a.g.l. are compared with two- 
hour average BND surface aerosol data centered on the flyby 
time. This comparison ties the airborne measurements to the 
surface measurements. Figure 6a and b show comparisons of 
AAO and BND submicrometer (aerodynamic diameter) and 
“total” size fraction aerosol light scattering coefficient data, 
respectively. In Fig. 6a, since the AAO submicrometer scat- 
tering data utilized the 1-1 nephelometer data at 545 nm, the 
BND submicrometer data from the 3-1 nephelometer were 
adjusted from 550 to 545 nm using the 450nm/550nm scat- 
tering Angstrom exponent. The plot shows that the measure- 
ments for Dp < 1 pm particles are quite similar (linear regres- 
sion slope forced through the origin is ~ 0.97) for surface 
and flyby data. The green dashed lines show similar popu- 
lation median values of 24.9 and 26.6 Mm -1 for the AAO 
and BND measurements, respectively. This suggests that for 
the atmospheric layer between the surface and ~ 240 m a.g.l., 
the aerosol loading of submicrometer particles was similar 
at least on a statistical basis and that the aircraft inlet ef- 
ficiently transmitted these particles to the instruments. Fig- 
ure 6b shows for the “total” aerosol population a regres- 
sion slope forced through the origin has a value of ~0.87. 
This suggests that some of the larger particles may not have 


been sampled by the aircraft instruments and/or that there are 
real differences in the aerosol loading of larger particles be- 
tween the surface and ~ 240 m a.g.l. The median values for 
the “total” size fraction comparison are similar at 29.9 Mm -1 
(AAO) and 31.1 Mm -1 (BND). Figure 6c and d show com- 
parisons of intensive parameters (the 550nm/700nm a s and 
550 nm a> o, respectively, for the “total” aerosol) measured 
on both platforms. In Fig. 6c, the median values of 2.04 
(AAO) and 2.18 (BND) agree well, and the linear regres- 
sion forced through the origin (slope of ~ 0.92) also suggests 
good agreement between the two sets of measurements. Fig- 
ure 6d shows a noisy relationship between the coq measure- 
ments at the surface and the flyby altitude, and this is pri- 
marily due to larger relative differences in the low er ap val- 
ues. The median m o values of 0.925 and 0.940 agree well 
however, and the slope of 0.988 suggests that these two sets 
of measurements are quite comparable with no major differ- 
ences observed between surface and the bottom of the lower 
column. The a s data indicate that for fine particles the surface 
and profile data sets agree well while the total size fraction 
data suggest either a very minor problem with the sampling 
of larger (i.e., supermicrometer) particles or real aerosol par- 
ticle size differences between the surface and the flyby alti- 
tude. While this exercise rules out major sampling problems 
with the airborne data, determining the potential contribution 
of each of these factors to the observed comparisons is be- 
yond the scope of this study. Esteve et al. (2012) performed 
a detailed analysis of the differences observed between the 
AAO in situ and column (sunphotometer) measurements and 
have eliminated some of the potential sources of discrepancy 
from consideration as being negligible. 

In order to relate the surface data to the AAO profile data 
we also compared the lowest flight segments conducted over 
the BND site with the lowest Lodge profile flight segments. 
Comparable flight segments from these two data sets were 
conducted within 15 min at locations centered ~ 15 km away 
from each other. The same parameters presented in Fig. 6 are 
compared again in Fig. 7. Excellent agreement is obtained for 
all extensive and intensive aerosol parameters at this altitude 
between the two locations; we therefore conclude that the 
AAO measurements made in the lowest profile flight segment 
are indeed representative of both the amount and the nature 
of the aerosol measured at and just above the BND site. 

3.1 Statistical distributions of aerosol properties aloft at 
low RH 

The statistical distributions of all AAO segment-average pro- 
file data (all LODG, TERA, AQB1, AQC1, and AQD1 pro- 
files conducted over the duration of the project) are shown 
as box-whiskers plots in Fig. 8. In these plots the median of 
each distribution is shown as a vertical line through the box, 
the ends of the boxes represent the 25th and 75th percentiles 
of the distributions, and the ends of the whiskers show the 5th 
and 95th percentiles. The data shown in Fig. 8 and the other 
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Fig. 6. Comparisons of BND surface aerosol data with concurrent AAO flyby data. Data points represent 5-min AAO low-level flight segment 
averages over the BND site vs. two-hour BND surface data centered on the flyby time, (a) <r S p for Dp < 1 pm particles, (b) cr sp for “total” 
size cut particles, (c) a s for “total” size cut particles, (d) cdq for “total” size cut particles. Red dashed line =1:1, green dashed lines show the 
median values for each population, and blue and black solid lines show the linear regression fit lines (black forced through the origin). 


multi-panel figures are from 399 of the 401 research profiles 
- two flights had power failures prior to reaching the profile 
location and no profile data were recorded. This represents a 
retrieval rate of > 99.5 % for the in situ aerosol data from the 
AAO program. The segment-average data in these plots have 
been adjusted from instmment conditions to ambient condi- 
tions of temperature and pressure, with the exception of the 
[CN], which are made using a volume flow rate measured 
at instrument temperature and pressure. All data reported in 
Fig. 8 are at instrument conditions of low (< 40 %) RH. 

For clarity, only the green wavelength in the multiple- 
wavelength instruments was used for presentation of the sta- 
tistical distribution plots in Fig. 8, but the spectral depen- 
dence of the measurements is shown by the a s and a a . Un- 
less otherwise noted, all plotted parameters were calculated 
from measurements either made at or adjusted to a common 
wavelength of 550 nm (the green wavelength of the TSI 3-A. 
nephelometer). The PSAP operating wavelength of 530 nm 
was adjusted to 550 nm by log-interpolation between the 530 
and 660 nm wavelengths of the PSAP for calculation of ex- 
tinction, single-scattering albedo, and radiative forcing effi- 
ciency. The wavelength of the 1 -/, nephelometer (545 nm) 
was adjusted to 550 nm using the segment-average a s from 
the TSI nephelometer (450 nm/550 nm) for calculation of the 
R s p . The a s data presented in Fig. 8 were calculated using 


the 550 and 700 nm wavelength scattering data from the TSI 
nephelometer, and the a a data were obtained using all three 
wavelengths from the PSAP. 

Figure 8 shows a clear decrease with altitude of all exten- 
sive aerosol parameters (scattering, absorption, extinction, 
[CN]). The medians of these extensive properties decrease 
more slowly with altitude above the 1830 m level, suggest- 
ing that the statistical top of the mixed layer for all flights 
during all seasons is probably near the 1830m a.s.l. level. It 
should be noted again that these flight levels are in fact pres- 
sure altitudes and therefore are approximate. 

The intensive aerosol properties in Fig. 8 are shown in 
panels g-1. Because the intensive properties are often ra- 
tios of the extensive properties, they can become quite noisy 
when aerosol extinction is low. For this reason, in this 
and all other multi-panel figures the intensive properties 
were limited to those flight segments that exceeded a min- 
imum threshold value of light scattering (segment-average 
Osp,550nm — 3.0Mm _l under dry (<40%RH) conditions). 
Segments with lower er sp , 550 nm values were excluded from 
the plotted intensive statistics presented in Figs. 8 through 
12. For lower altitude segments (i.e., those conducted at or 
below the 1520m a.s.l. level in the boundary layer), only a 
few segment-averages from the 399 valid profiles fell below 
the light scattering threshold. For the higher-altitude flight 
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Fig. 7. Comparisons of AAO lowest-segment profile data with BND tower flyby data. Data points represent 5-min flight segment averages of 
aerosols at the same altitude and displaced by ~ 15 km. (a) cr sp for D p < 1 pm particles, (b) cr sp for “total” size cut particles, (c) a s for “total” 
size cut particles, (d) coq for “total” size cut particles. Red dashed line =1:1, green dashed lines show the median values for each population, 
and black solid line shows the linear regression fit line (forced through the origin). 


levels, the exclusion rate ranged from 23 % at 1 .8 km to 59 % 
at 4.6 km of the valid flight segments. Our free troposphere 
(FT) aerosol statistics are therefore biased toward the in- 
tensive aerosol properties present during less pristine condi- 
tions, and it is unknown how these would compare with the 
same properties derived during clean conditions. 

The intensive properties generally show less variability in 
the vertical than do the extensive properties, which is broadly 
consistent with previous studies (e.g., Sheridan and Ogren, 
1999; Sheridan et al., 2002; Andrews et al., 2004, 2011). 
Median coq (Fig. 8g) varies from ~ 0.925 at the lowest flight 
level to ~ 0.94 at the statistical top of the mixed layer, and 
stays relatively constant above that. The differences between 
the median a > o measurements at low and high altitudes in this 
study (~ 0.02, or roughly a few percent relative difference) 
are within the realm of measurement uncertainty given the 
estimated total measurement uncertainties in the nephelome- 
ter (~ 5-10 %) and PSAP (~ 10-20 %). Scattering Angstrom 
exponent (Fig. 8j) decreases with altitude above the bound- 
ary layer, an observation consistent with the mid-Atlantic 
coast study of Taubman et al. (2006). Median top of atmo- 
sphere (TOA) radiative forcing efficiency (RFE, Sheridan 
and Ogren, 1999) varies from —25 to — 28 W m' 2 per unit 


AOD throughout the column, indicating that, at least in a sta- 
tistical sense, the surface measurements can in some cases be 
used to estimate direct aerosol radiative forcing if a suitable 
aerosol optical depth measurement is available. 

The spectral dependence of aerosol light absorption in 
the atmosphere depends upon the absorption spectra of the 
individual aerosol constituents and also their size, shape 
and/or mixing state (Bergstrom et al., 2007). A larger a a 
suggests the possibility that organic compounds may be on 
or in the particles to a larger extent (Kirchstetter et al., 
2004; Andreae and Gelenscer, 2006; Bergstrom et al., 2007). 
Schnaiter et al. (2006) showed that by altering the burn con- 
ditions and producing propane soot particles with varying or- 
ganic/total carbon (OC/TC) ratios, the a a varied from near 1 .0 
at ~ 8 % OC/TC to ~ 3.5 at ~ 50 % OC/TC. The median a a 
(Fig. 8k) from the mixed atmospheric aerosols sampled dur- 
ing the AAO program increased from just above 1.0 in the 
lower portion of the mixed layer to ~ 1.4 at the highest flight 
altitude. 

Observations of elevated a a values are not, however, lim- 
ited to organic aerosols. Bergstrom et al. (2004), Collaud 
Coen et al. (2004), and Fialho et al. (2005) discussed the 
wavelength dependence of absorption for desert dust or 
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Fig. 8. Cumulative statistics of measured and calculated parameters for all AAO research profiles. Aerosol measurements are for dry 
(RH < 40 %) conditions and are adjusted to ambient temperature and pressure, with the exception of the CN concentration, which is re- 
ported at instrument temperature and pressure (volume flow rate used). 


dust/pollution aerosol mixtures. These studies found elevated 
a a values in the range of 1.5 to nearly 3.0 for aerosol with a 
significant dust component. We cannot distinguish from the 
optical measurements whether the observed a a values at high 
altitude are from dust or organic aerosols, or a mixture of 
both. 

The fact that portions of the distributions at all altitudes 
show a a values less than 1 is interesting given that the theo- 
retical small particle limit for light absorption (spherical par- 
ticles, constant index of refraction) has been shown to have 
a A. -1 dependence (Bergstrom, 1973; Bohren and Huffman, 
1983). This could have to do with measurement uncertain- 
ties associated with the small light absorption coefficients 
(the threshold for inclusion of segment averages in the statis- 
tics was er S p,550nm > 3. 0 Mm~ 1 , which generally resulted in a 
cLp.ssonm threshold of around 0.2 Mm~ 1 - this threshold may 


be too low for practical calculation of a a ). It could however 
be partially explained by the fact that the imaginary part of 
the index of refraction increases slightly with wavelength for 
some types of particles and this could cause the a a to dip be- 
low the theoretical limit of 1 .0 (Bond and Bergstrom, 2006, 
and references therein). 

The profiles of R S p, b, and a s all suggest that there is 
a larger relative contribution to scattering extinction from 
larger particles with increasing height above the surface up 
to 4.6 km. This could be explained by more sources of small 
particles (e.g., combustion aerosols, photochemical smog) 
at the surface, long-range transport of supermicrometer par- 
ticles aloft, or selective wet scavenging processes at al- 
titude. Less-hygroscopic aerosols may have different size 
distributions than water-absorbing aerosols or may contain a 
larger fraction of OC, and organic aerosols have been shown 
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Fig. 8. Continued. 

to resist wet removal processes (McFiggans et al., 2006; and 
references therein). Similarly, uncoated dust or dust coated 
with organics would presumably not be very hygroscopic and 
could remain aloft for extended periods. 

The in situ aircraft data for the statistical distributions of 
aerosol properties presented in Fig. 8 have been plotted in 
time-altitude space in Fig. 9. These contour plots show how 
the aerosol properties vary with altitude over the duration of 
the AAO program. In these plots, individual flight segments 
are represented as black dots, and the segment-average val- 
ues for each parameter have been contoured. The variation 
in sampling altitudes is caused by a seasonal variation in the 
relationship between pressure and geometric altitudes. The 
extensive parameters (Fig. 9a-d) show a wavelike pattern, 
with relatively larger values of these parameters extending 
to higher altitudes in the summertime. This may be because 
of the increased height of the top of the mixed layer during 
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the warmer months due to enhanced vertical mixing, signifi- 
cant seasonal increases in the aerosol sources (primary and/or 
secondary) beginning in the spring and lasting through fall, 
and/or fewer or less efficient removal mechanisms at work 
during this period. 

The annual ambient temperature and RH cycles are 
shown in Fig. 9e and f, respectively. A ~4 month gap 
in the ambient RH, from late 2006 through early 2007 
was caused by a malfunctioning RH sensor. As in Fig. 8, 
the intensive parameters have been plotted in panels 9g-l, 
with the same noise-filtering threshold of segment-average 
°sp,550nm > 3.0Mm -1 . Panel 9g shows a weak annual cycle 
in a>Q, with the lowest values in the boundary layer in the fall. 
This supports the September-October minimum in the long- 
term a>o record observed at the Bondville surface monitoring 
site (Delene and Ogren, 2002), as well as the similar fall co o 
minimum at the mid-continental SGP site (Sheridan et al., 
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Fig. 9. Time series contour plots of AAO measurements. Black dots represent individual flight segments. 


2001). Consistent with Fig. 8, the data displayed in panels 
9h-j suggest that larger particles have a greater influence on 
the aerosol optical properties at higher altitudes. 

The contour plots of intensive properties show an anoma- 
lous high altitude feature in early 2009 that is visible in the 
upper right portion of most of the panels. Panels 9h-j sug- 
gest that larger particles were present at this time, and too 
was higher than typical, with values exceeding 0.94. Over- 
all RFE was lower than that typically observed for this pro- 
gram at < —25 W m -2 . Larger, mostly scattering aerosols at 
higher altitudes could be long-range intercontinental trans- 


port of Asian dust and/or pollution or aged biomass smoke 
advecting over the central US during the springtime (Van- 
Curen, 2003; Augustine et al., 2008). The springtime Asian 
aerosol, which has both a dust and pollution component, has 
been shown to stay aloft for very long periods of time and to 
travel great distances (Rahn et al., 1977). 

Median profiles of aerosol properties segregated by sea- 
son are shown in Fig. 10. Seasons are described here accord- 
ing to the conventional three-month definitions of spring = 
MAM, summer = JJA, fall = SON, and winter = DJF. In all 
plots, colored lines show the spring, summer, fall, and winter 
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Fig. 9. Continued. 

median profiles along with the black annual median profile. 
Black horizontal lines show the extent of the 25th and 75th 
percentiles of the annual distributions. The seasonal plots for 
the extensive properties confirm that for most of the sam- 
pled vertical column, the largest aerosol loadings over the 
BND site were in the summer season and the lowest were 
during the winter. Above 3kma.s.l„ the median spring ex- 
tinction values exceeded those for all other seasons. This sug- 
gests the occasional presence of elevated aerosol layers in the 
springtime over the BND site. When interpreted in conjunc- 
tion with the contour plots, these data suggest that springtime 
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aerosol incursions at high altitude over central Illinois are in- 
frequent and episodic in nature rather than a regular event 
each spring. 

The seasonal plot of cdq supports the contour data pre- 
sented in Fig. 9. The lowest &>o values at low altitude ap- 
pear in the fall, and the seasonal values generally fall be- 
tween ~ 0.92-0.95. The wintertime 3.0 and 3.7 km levels 
show slightly larger coq median values of ~ 0.96, but this may 
be an artifact of the limited sample population. Fewer than 
25 % of the 65 wintertime flight segments at these high al- 
titudes exceeded the er S p, 55011 m threshold value. Panels lOh-j 
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Fig. 10. Median seasonal profile plots of AAO measurements. Black profile in each plot represents all aircraft profiles. Horizontal bars show 
the extent of the 25th and 75th percentiles of the total distribution. 


show the general trends of decreasing A* S p, b, and a s with in- 
creasing altitude, consistent with the presence of larger par- 
ticles at higher altitudes. 

Cumulative seasonal plots of the median scattering, ab- 
sorption, and extinction (all at 550 nm), and [CN] are pre- 
sented in Fig. 11. These plots of the extensive parameters 
show the percentage of the total column amount of each pa- 
rameter that occurs below a given altitude in the atmospheric 
column from the surface to 4.6 km. For example, the cumu- 
lative extinction plot shows that for all seasons ~ 80 % of the 
extinction in the column is observed at or below 1.8 km al- 
titude, while in the winter ~ 90 % of the column extinction 
is below that altitude. In general, winter observations show 
that a larger fraction of the total column aerosol is present at 
lower altitudes than during the other seasons. This is prob- 
ably because of the lower mixed layer heights in the winter 


over this site. The green lines representing spring medians 
lie above the other seasonal lines for scattering, absorption 
and extinction, and can probably be explained by the fact 
that high altitude (3 km altitude and above) aerosol layers 
were present more frequently during the springtime, so that a 
smaller percentage of the total column abundance of these 
parameters exists at the lower altitudes. These higher me- 
dian values can be seen clearly in the seasonal profile plots 
in Fig. 10. 

3.2 Aerosol scattering hygroscopic growth 
measurements 

Over the first two years, we had difficulty for a variety of rea- 
sons in running the humidified nephelometer system. These 
problems included intermittent electrical power and signal 
connections, drifting calibrations, overnight cold-soaking of 
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Fig. 10. Continued. 

the humidifiers and nephelometers (leading to water freez- 
ing in the system and ruptured humidifier membranes), and 
an inability to accurately control the internal RH of the hu- 
midified nephelometers at their respective set points which 
led to condensation of water vapor inside the instruments. 
Our program was designed to have a pilot fly an automated 
airborne aerosol system and do very limited pre-flight work, 
so many of these problems, including the drifting calibra- 
tions, were not fixable without service visits. The problem of 
proper control of the nephelometer internal RH was finally 
solved late in the program by our development of custom, 
real-time software PID control. Hygroscopic growth data for 
the early flights is not retrievable. 

Valid aerosol hygroscopic growth data were first obtained 
in the fall of 2008. During this period the hardware issues 
were gone but the problem of drifting calibrations for the hu- 
midified nephelometers remained. We visited the AAO sev- 




Scattering Angstrom exponent 



Radiative forcing efficiency (W m ! AOD ') 


eral times over this year and recalibrated the instruments, and 
typically flew one to several flights before the calibrations 
drifted to outside acceptable limits. For unknown reasons, 
these calibration drifts were not gradual but relatively large 
step changes between flights. Over the last year of the project 
we obtained /(RH) data from 12 profiles where calibra- 
tions were good (based on frequent calibration checks) and 
all nephelometers were functioning properly. These data are 
shown in Fig. 12. As with the other intensive measurements, 
the reported hygroscopic scattering growth data were lim- 
ited to segments where the segment-average, dry <r S p, 550 nm 
was > 3.0 Mm~'; this limited the data to 80 total flight seg- 
ments at the various altitudes. This constraint significantly 
reduced the number of valid high-altitude segments we report 
but should reduce the statistical uncertainty inherent in ratios 
of very small numbers. Figure 12 illustrates that the distri- 
butions of /(RH) and y in the FT are clearly different from 
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Fig. 11. Seasonal median plots showing the cumulative percentage of each extensive parameter in the lower atmospheric column (surface to 
4.6 km a.s.l.) that occurs at or below a given altitude. 
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Fig. 12. Cumulative profile plot of the statistical distributions of 
/(RH). The number of flight segments at each altitude that com- 
prise the distributions is shown at right. The gamma values were 
derived by fitting curves as described in Table 2 to the humidified 
nephelometer measurements. The /(RH) is the scattering hygro- 
scopic growth factor from 40 to 85 % RH. 


those in the planetary boundary layer (PBL, < ~ 1.8 km a.s.l. 
in the profiles shown in Figs. 8 and 10), with the lower val- 
ues (/(RH) < 2.0) observed in the FT. This could be because 


of surface sources of hygroscopic particles and/or because of 
the role cloud-processing plays in removing a fraction of the 
more hygroscopic material at the higher altitudes, leaving be- 
hind an aerosol enriched in less hygroscopic particles (e.g., 
Weingartner et al., 1999; Marcq et al., 2010; Berkowitz et al., 
2011). 

Esteve et al. (2012) have already presented a detailed 
discussion on the limited AAO /(RH) measurements and 
their comparison with previous hygroscopic growth studies 
at BND so only a brief synopsis is presented here. Aerosol 
chemistry data from selected periods during the Koloutsou- 
Vakakis et al. (2001) study period yielded a y value of 0.67 
by utilizing the Quinn et al. (2005) parameterization for 
/(RH), which is based on the relative amounts of particu- 
late organic matter and sulfate in the aerosol. A y of 0.67 re- 
sults in an estimated /(RH) (at 82.5 % RH) of 2.28 using the 
Quinn et al. (2005) parameterization, which is higher than the 
/(RH) value of 1.5 (effective y = 0.33) derived from the hu- 
midified nephelometry measurements of Koloutsou-Vakakis 
et al. (2001). Note that the Koloutsou-Vakakis /(RH) mea- 
surements were based on a high-RH scattering measurement 
at 82.5 %, which is slightly lower than the 85 % employed in 
this study. Independent aerosol chemistry data from the IM- 
PROVE (Interagency Monitoring of PROtected Visual Envi- 
ronments) network site at BND between January 2006 and 
December 2008 also gives a y value of 0.67 (/(RH) = 2.28) 
using the Quinn et al. (2005) parameterization. A median y 
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Fig. 13. Cumulative box/whiskers plot of AAO aerosol light scatter- 
ing profile data adjusted to ambient RH. Individual segment average 
scattering data were adjusted using the median y value at each alti- 
tude from Fig. 12 and the segment average RH. Red line represents 
the median dry values from Fig. 8a. 


value of 0.51 for the AAO aircraft measurements was de- 
rived by fitting curves as described in Table 2 to the humid- 
ified nephelometer measurements from the 80 valid AAO 
flight segments where all humidified nephelometers were 
working properly. This value of y results in an estimated 
/(RH) (at 85 % RH) of 2.03, intermediate to the values de- 
rived from the chemical composition and Koloutsou-Vakakis 
approaches (Esteve et al., 2012). For comparison, the me- 
dian boundary layer (all flight segments conducted below 
~ 1.8 km altitude) y value was 0.53 (/(RH) = 2.08) and the 
median free troposphere (all flight segments above ~ 1.8 km 
altitude) y value was 0.44 (/(RH) = 1.84). 

The cumulative dry aerosol <7 sp distributions shown in 
Fig. 8 were adjusted to ambient RH and are presented 
in Fig. 13. Individual flight segment-average cr sp values at 
550 nm were adjusted to ambient RH using the median 
altitude-specific y value for all valid segments and the 
segment-average RH. The median y value for each flight 
level was used because so few flight segments were con- 
ducted with valid hygroscopic growth measurements. It is 
implicit in this analysis that the median y values at each 
flight altitude adequately represent the hygroscopic growth 
characteristics of the aerosols during flights when the hu- 
midified nephelometry system was not working properly. We 
observed no obvious long-term changes in aerosol character- 
istics over the duration of the program to lead us to believe 
that this assumption was not valid. Additionally, the segment- 
average RH was typically low enough that small to moderate 
changes in the y parameter would not result in major changes 
in the adjusted er sp values. 

Even though the adjustment of scattering to ambient RH 
is nonlinear in RH, we felt using the segment average rather 
than the instantaneous RH was appropriate for two reasons. 


First, the RH measured on most of our flight segments was 
not very high, and RH variability at lower RH does not cause 
large changes in the optical properties. Figure 8f shows that 
median RHs for flight segments in the mixed layer were 50- 
60% and were lower at higher altitudes. Only a small frac- 
tion of all flight levels show segment-average RH values at 
or above 70 %, where the growth function becomes steeper. 
Secondly, RH variability was consistently small over the 5 
and 10 min horizontal segments flown over the course of the 
program. Considering all AAO flight segments, the average 
standard deviation (SD) of RH measurements during indi- 
vidual flight segments was 2% (with a SD of the SD val- 
ues of ±2 %). That is, on average, each AAO flight segment 
showed a relatively small variability (expressed as the SD) of 
2 % RH. While there were a few flight segments conducted at 
higher RH that also showed a larger RH variability, this was 
not a commonly observed feature of the data set. 

In general, the changes from the dry scenario for light scat- 
tering shown in Fig. 13 are relatively small, with median am- 
bient scattering values in the PBF higher than the dry case 
(dry median values shown by red line) by only a few Mm - 1 . 
This can be explained by the fact that flights typically oc- 
curred during sunny to partly cloudy conditions where am- 
bient RH was relatively low (typically <60%). Figure lOf 
shows that even during the summer, median RH’s in the PBF 
during AAO flights were ~ 57-65 %. Relative humidities in 
this range are not high enough to dramatically increase the 
dry scattering values. No attempt was made to scale AAO 
light absorption data (or single-scattering albedo and extinc- 
tion data) due to a lack of information on the hygroscopic 
growth function for er ap . The in situ AAO extinction mea- 
surements used for comparison with the satellite-derived ex- 
tinction reported below have thus had their scattering com- 
ponent adjusted to ambient RH conditions while the absorp- 
tion component was held at low RH conditions. These in situ 
measurements therefore most likely represent a lower limit 
to the true er e p values at ambient RH, assuming a ap increases 
with increasing RH. Since RH was not high in most flight 
segments and since <r ap is usually a small fraction of cr ep , the 
reported AAO <r ep values should be close to the true ambient 
values. 

3.3 Satellite vs. in situ extinction comparisons 

One of the objectives of this program was to obtain in situ 
aerosol data during satellite overflights for potential satel- 
lite retrieval verification. The AAO <r ep data used in all com- 
parisons with the CAFIPSO lidar have been adjusted as dis- 
cussed above to ambient conditions of T, p, and RH to the 
extent possible for these comparisons. The scattering com- 
ponent of AAO extinction for each comparison was adjusted 
to ambient RH utilizing the median y value for that specific 
altitude, and was then adjusted to the 532 nm wavelength 
of the CAFIPSO lidar based on the 450/550 nm a s . The ab- 
sorption component of AAO extinction was not adjusted to 
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ambient RH given the lack of understanding of the RH- 
dependence of er a p in ambient aerosols. The in situ absorp- 
tion measurement was adjusted to the lidar wavelength using 
the 530nm/660nm a a . 

The three primary A AO profile locations were chosen to 
be within 10 km of three A-Train satellite tracks (Fig. 3), and 
as shown in Fig. 7 and reported in Anderson et al. (2003), 
the aerosol properties are expected to be similar over this 
distance. Over the course of the program 63 AAO flights 
were conducted at these locations where the aircraft was in 
the PBL within two hours of the CALIPSO satellite over- 
pass. It is reasonable to ask whether there were any system- 
atic differences between the aerosol properties measured dur- 
ing the satellite and non-satellite flights. One could surmise 
that, since the satellite underflights were constrained opera- 
tionally to occur on less cloudy days, the aerosol properties 
could be different in drier conditions. It was observed, how- 
ever, that there were not major statistical differences in either 
the extensive or intensive aerosol properties between flights 
with and without satellite overpasses. For example, the me- 
dian AAO-measured a ep and coo values (at 550 nm) measured 
during the 63 overpass flights at the five segment altitudes in 
the boundary layer agreed to within 2 Mm -1 and 0.005, re- 
spectively, with the analogous measurements from the 336 
non-overpass flights. 

CALIPSO Level 2 aerosol extinction profiles are retrieved 
based on 5, 20, and 80km horizontal averages of total at- 
tenuated backscatter and are reported every 5 km along the 
track at 60 m vertical resolution. The comparisons with in 
situ aerosol extinction are based on the AAO measurement 
being matched with the CALIPSO retrieval from the closest 
60-m range bin to the mean flight segment altitude. For this 
comparison, 35-km (i.e., 7 extinction profiles) along-track 
averages centered on the closest profile were used so that the 
satellite would sample a similar airspace as the airplane. We 
initially analyzed the CALIPSO data both using 180-m ver- 
tical averaging and also with no averaging, and the profile 
statistics did not change very much. We decided not to av- 
erage vertically to avoid incorporating aerosol not observed 
by the aircraft into CALIPSO profile data, but acknowledge 
that layer heights may change by 120-180 m over a distance 
of 35 km. As discussed below, this is not believed to be a 
prevalent feature in this data set. 

CALIOP lidar aerosol extinction data (Level 2 Aerosol 
Profile Product, Version 3.01 processing) were retrieved for 
28 of the 63 collocated profiles. CALIPSO data are quality- 
screened based on metrics used to generate CALIPSO Level 
3 aerosol profile products from Level 2 aerosol extinction 
data (Winker et al., 2012) and are described in Appendix A. 
Quality screening removed approximately 4 % of all possi- 
ble collocated CALIPSO extinction samples. “No retrieval” 
cases most likely mean that it was either too cloudy, there 
was not enough aerosol present, and/or there was too much 
noise in the daytime CALIPSO profile to retrieve extinc- 
tion at the resolution required. Inspection of the CALIPSO 


browse images (using the vertical feature mask) revealed that 
for some of these 35 “no retrieval” cases, CALIPSO de- 
tected aerosols farther away from the aircraft profile loca- 
tion but not within the 35-km comparison window even for 
clear sky conditions. For other flights, overlying clouds atten- 
uated the lidar and obscured the viewing beneath, so it is not 
clear whether detectable levels of aerosols were present in 
the boundary layer in these cases. The 28 AAO flights where 
CALIPSO also retrieved extinction showed in situ median 
aerosol extinction up to about 20 % higher than the 35 flights 
where CALIPSO observed no aerosol. Medians of the distri- 
butions of 550 nm extinction for the 35 “no retrieval” cases 
were typically in the 15-35 Mm -1 range for flight levels in 
the boundary layer, and as discussed below this suggests that 
CALIPSO would not have retrieved aerosol extinctions in 
this range a substantial fraction of the time. We limited our 
retrievals to daytime profiles because the nighttime lidar pro- 
files, while considerably less noisy, were too far removed in 
time for valid comparisons. For example, a 12-h time differ- 
ence with winds at 10 km per hour puts the aerosols 120 km 
off track at AAO flight time, which is too far away to expect 
the PBL aerosols would be similar (Anderson et al., 2003). 
Additionally, diurnal changes in the height of the boundary 
layer would also affect the shape of the vertical profile. 

Figure 14 shows statistical comparisons of AAO in situ 
extinction and CALIPSO lidar extinction data for the 28 pro- 
files where collocated comparisons were possible. In order 
to understand the effects of boundary layer clouds on the 
comparisons, separate analyses of the cases with and with- 
out boundary layer clouds were performed. Of the 28 flights 
where CALIPSO detected aerosols in the proximal 35-km 
satellite track, the vertical feature mask images showed 8 
cases without and 20 cases with boundary layer clouds. The 
mean extinction profiles for the 8 cases without boundary 
layer clouds (Fig. 14a), the 20 cases with boundary layer 
clouds (Fig. 14b), and all 28 cases (Fig. 14c) are presented. 
The shaded envelopes represent standard deviations of the 
observations and the CALIPSO error bars are the propagated 
uncertainties of the lidar measurement (Winker et al., 2009; 
Young and Vaughan, 2009). The mean extinction for the 
cases without boundary layer clouds is smaller below ~ 3 km 
for both instruments compared to the cases with boundary 
layer clouds. The agreement below 1 km is better between 
AAO and CALIPSO for the case without boundary layer 
clouds, except for the spike in CALIPSO extinction popu- 
lation at 0.9 km. The extinction browse images show gran- 
ules on two flights that could have caused this spike in the 
mean extinction profile. We suspect unidentified boundary 
layer clouds given the altitude, although this will require 
more study. It is interesting that the spike at 1.2 km in the 
case with boundary layer clouds is reduced significantly in 
the case without boundary layer clouds. This supports our 
hypothesis that cloud contamination could be occurring, par- 
ticularly at that altitude. 
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Fig. 14. Cumulative statistics of the mean extinction profiles of the 28 collocated in situ AAO and CALIPSO lidar comparisons. The shaded 
envelopes are standard deviations of the observations and the CALIPSO error bars are the propagated uncertainties of the lidar measurement, 
(a) No observed boundary layer clouds in proximal satellite track (8 cases), (b) Boundary layer clouds observed in proximal satellite track 
(20 cases), (c) All 28 collocated flights where CALIPSO retrieved extinction in the profile. 
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Fig. 15. Individual flight segment extinction comparisons, color coded by (a) altitude and (b) CALIOP lidar ratio. 


These statistical compilation plots suggest that CALIPSO 
tends to overestimate the extinction in the boundary layer at 
this location while underestimating the extinction in the free 
troposphere. Imperfect cloud masking may explain some of 
the higher satellite extinction observations that influence the 
larger mean values between ~ 0.4 and 1 .5 km altitude. Some 
individual lidar profiles show large “extinction” peaks that 
were observed to be much smaller by aircraft at altitudes 
where scattered clouds were present. The lower CALIPSO 
extinction values in the FT are likely related to CALIOP 
aerosol detection limits. In the flight segments where AAO 
detected relatively low levels of in situ extinction, CALIOP 
generally did not get a valid retrieval (defined as “clear air” 
cases). In these cases extinction values of 0.0 Mm' 1 were av- 
eraged into the compilations, so these free tropospheric lidar 
values are clearly lower limits. For example, the CALIPSO 
mean value at 3.7 km is 0.0 Mm -1 , meaning that CALIPSO 
did not get a valid retrieval at this altitude for any of the 


28 collocated cases. This is not surprising given the low in 
situ extinction values and the noisy daytime retrievals. There 
is currently no known uncertainty to assign for CALIPSO- 
defined clear air cases. Limiting the averages to the aerosol 
layers detected by CALIPSO would mean that the averages 
only represent enhanced conditions and would be biased 
high. 

Figure 15 shows the extinction comparison for each flight 
segment color coded by altitude (a) and lidar ratio (b). In 
Fig. 15a, the data points lying along the x-axis (CALIPSO 
extinction = 0.0 Mm -1 ) are mostly from the higher-altitude 
comparisons where in situ extinction was low and aerosol 
layers were undetected by CALIPSO. There are a few data 
points that show higher in situ extinction that were not de- 
tected by the CALIOP lidar. There is also a population of 
large (> 100 Mm -1 ) CALIPSO extinction values, typically 
at lower altitudes, that do not correspond with equally large 
AAO extinction values. These most likely represent the 
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Table 4. Aerosol types and lidar ratios identified by CALIPSO dur- 
ing AAO overpasses. 


Aerosol Type 

Lidar Ratio 
(sr) 

Percent of 
Samples (%) 

Dust 

40 

15 

Polluted dust 

55 

46 

Polluted continental 

70 

24 

Smoke 

70 

15 


comparisons where the cloud screening process should have 
eliminated these satellite data from consideration but did not. 
This population of large CALIPSO extinction values rela- 
tive to AAO could, however, result from either misidentifi- 
cation of aerosol type or incorrect assignment of lidar ratio, 
or both. The lidar ratio (ratio of extinction to backscatter) 
is required by the CALIPSO algorithms to retrieve extinc- 
tion from backscatter and varies by aerosol type. Figure 15b 
shows the AAO vs. CALIPSO extinction with the lidar ratio 
identified for each comparison point. Many of these higher 
CALIPSO extinction values were derived using high lidar ra- 
tios of around 70 (dark red circles). 

Aerosol types identified by CALIPSO are shown in Ta- 
ble 4 for the 28 collocated flights. Most aerosol samples in 
this regime (46 %) are classified as polluted dust with a li- 
dar ratio of 55 sr. However, 39 % of samples are classified 
as either polluted continental or smoke aerosol types, both 
having lidar ratios of 70 sr. If this lidar ratio is too large for 
aerosol in these cases the retrieved extinction will be overes- 
timated and could contribute to the large CALIPSO extinc- 
tion values since errors in retrieved extinction are roughly 
proportional to errors in the assumed lidar ratio. Indeed, An- 
derson et al. (2000) measured lidar ratio values at BND of 
64 ± 4 sr which, if representative of aerosol measured during 
these 28 flights, could account for some of the discrepancy 
in extinction below 1.5 km. Aerosols classified by CALIPSO 
as dust comprise 15 % of the samples and the vast majority 
of higher altitude (> 3.0km) samples were classified as this 
aerosol type. This supports the suggestion from the in situ 
AAO intensive parameters of larger mean particle sizes at 
higher altitudes. 

Our long-term aircraft data permit the compilation of 
CALIPSO extinction detection efficiency as a function of in 
situ extinction level. The data from all 28 collocated aircraft 
profiles are shown in Fig. 16. The number of discrete com- 
parison points represented in this figure is 244; that is the 
number of the 280 total flight segments that the AAO made 
a valid extinction (i.e., scattering plus absorption) measure- 
ment. The in situ AAO extinction is grouped into 10 Mm -1 
bins, and the number of cases in which CALIPSO detected 
extinction in each range bin is shown above each column. 

Figure 16 can provide a clue as to how prevalent the 
potential problem of aerosol layers changing altitudes over 



AAO Extinction (Mm -1 ) 


Fig. 16. CALIPSO extinction detection frequency as a function of 
AAO in situ extinction. Numbers above bars indicate the number 
of valid CALIPSO aerosol extinction retrievals falling within each 
AAO extinction range bin. 


the 35-km path might be. For in situ extinction values 
above 50 Mm -1 , CALIPSO retrieved a valid extinction value 
~ 95 % of the time. The only bin where CALIPSO clearly 
missed the extinction was in the 140-150 Mm -1 bin, where 
it retrieved only 2 of the 3 segments that fell into that 
bin. Closer inspection of this one missed case revealed that 
CALIPSO retrieved a clear air extinction value near what 
was measured by the aircraft in the adjacent vertical range 
bin. The CALIPSO measurement for this point was made 
22 min after the AAO measurement. It is possible that the 
aerosol layer rose 30-60 m over this period of time and was 
classified into the adjacent bin. There are other possibili- 
ties also, but a detailed analysis of these individual com- 
parisons is beyond the scope of this paper. In this one case, 
we can clearly say that vertical averaging of the CALIPSO 
data would have helped to detect a layer sampled by the air- 
craft. This does not appear to be a problem that was observed 
frequently. It is possible that CALIPSO also occasionally 
missed retrieving extinction in the 60-70, 70-80, and 80- 
90 Mm" 1 bins, but this could also be due to cloud contami- 
nation issues as much as aerosol layers changing height. 

For in situ extinction levels larger than 50 Mm" 1 , 
CALIPSO makes a retrieval ~ 95 % of the time (77 out of 
81 cases). For the 40-50 Mm -1 extinction bin, CALIPSO’s 
retrieval frequency drops to < 80 %. Below about 20 Mm" 1 , 
CALIPSO detected aerosol extinction only in about 1 1 % (12 
out of 109) of the cases, and below 10 Mm -1 , CALIPSO 
rarely (1 case in 76) retrieved any extinction. For all 244 
cases, a 50 % probability of detection falls at an in situ ex- 
tinction level of 20-25 Mm - 1 . 
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4 Conclusions 

Over 400 vertical profiles of in situ aerosol properties were 
obtained using a light aircraft over eastern and central Illi- 
nois between June 2006 and September 2009. A majority of 
the aircraft profiles were flown near a NOAA/ESRL surface 
aerosol monitoring station where identical aerosol measure- 
ments were being conducted. The vertical profile measure- 
ments were conducted to build a robust database of aerosol 
properties aloft for the purposes of assessing aerosol variabil- 
ity in the vertical and determining how representative sur- 
face aerosol measurements are of those in the lower column. 
In addition, individual profiles from this program were com- 
pared with CALIPSO lidar extinction measurements in order 
to provide “air truth” measurements for satellite retrieval val- 
idation exercises. 

Flyby comparison measurements showed good agreement 
between those from the low altitude aircraft flyby segment 
and those conducted with similar or identical instruments at 
the surface. Excellent agreement was obtained for light scat- 
tering by submicrometer aerosols, and agreement was ac- 
ceptable even for the “total” size fraction (particles up to the 

5 pm diameter efficiently sampled by the aircraft inlet), indi- 
cating no major sampling problems for the aircraft inlet or 
instruments. 

Extensive aerosol properties show that the amount of 
aerosol decreases strongly with increasing altitude to the top 
of the mixed layer, and more slowly after that. Much less 
variability was observed in the intensive properties. Several 
aerosol parameters suggest a lesser contribution to aerosol 
optical properties from smaller particles at the higher alti- 
tudes. This could be due to increased small particle sources 
at the surface or more efficient removal mechanisms for the 
smaller particles at higher altitudes owing to their size and/or 
composition. The median aerosol absorption Angstrom ex- 
ponent increased from just over 1.0 in the lower portion of 
the mixed layer to ~ 1 .4 at the highest flight altitude. A larger 
a a suggests that dust or organic material may be on or in the 
particles to a larger extent. A larger dust or organic compo- 
nent in the aerosol may be responsible for the larger parti- 
cles surviving precipitation and scavenging processes due to 
a lower hygroscopicity of those substances. Aerosol scatter- 
ing hygroscopic growth measurements made during a subset 
of profiles show less-hygroscopic aerosols at higher profile 
altitudes. Scattering hygroscopic growth factor medians were 
observed to be > 2.0 in the PBL and < 2.0 in the free tropo- 
sphere. 

Aerosols over Illinois show clear statistical differences 
among seasons. Seasonal differences include more extinc- 
tion in the column in the summer and less in the winter, 
with spring and fall profile medians falling in between the 
extremes. The springtime profiles, however, show the largest 
extinction values of any season above 3kma.s.l. In general 
the lowest o>q values in the column were observed in the fall, 
and this is consistent with the fall minimum in a>o observed in 


the surface measurements. AAO-derived RFE showed very 
little vertical or seasonal variability. Barring large differences 
in aerosol hygroscopicity and RH (which are observed in 
some individual profiles but not in the cumulative statistics), 
these data suggest that surface measurements of aerosol opti- 
cal properties could be used, on average, to estimate column 
RFE. However, an instantaneous retrieval of RFE using sur- 
face properties could be in error. 

The individual AAO/CALIPSO extinction comparison 
cases that comprised this study hold a wealth of informa- 
tion for verifying the satellite retrievals. Our purpose was 
to show how well, in a statistical sense, the satellite data 
agreed with our in situ extinction data. AAO underflights of 
the CALIPSO satellite show generally good agreement be- 
tween aircraft-measured and satellite-retrieved extinction for 
our limited subset of collocated profiles where clouds were 
not prevalent. We have observed that CALIPSO often overes- 
timates extinction in the PBL at this location, possibly due to 
interferences from clouds that are not screened properly from 
the retrieval, but determining the specific causes for these 
discrepancies in individual profiles is beyond the scope of 
this work. Our long-term data suggest that if light extinction 
at 532 nm in atmospheric aerosol layers exceeds 50 Mm -1 , 
CALIPSO will retrieve an extinction value ~ 95 % of the 
time. The in situ extinction threshold at 532 nm for a 50 % 
probability of detection by CALIPSO is 20-25 Mm~ 1 . We 
expect to investigate the individual profile comparison cases 
in more detail in a future paper. 


Appendix A 

CALIPSO data quality screening 

CALIPSO aerosol extinction data are quality screened us- 
ing metrics reported in the Level 2 aerosol profile data prod- 
uct which indicate the algorithm’s ability to determine the 
correct atmospheric feature type (aerosol or cloud) and the 
final state of the extinction retrieval algorithm. These qual- 
ity screening filters are based on those used to generate 
the CALIPSO Level 3 aerosol profile product from Level 
2 aerosol profile data; rationale and further details are de- 
scribed in Winker et al. (2012). 

Each layer detected by the CALIPSO feature classification 
algorithms is given a Cloud-Aerosol-Discrimination (CAD) 
Score, a metric indicating the confidence in determining if 
the layer is an aerosol or a cloud. If CAD = — 100, it reflects 
that layer is most likely aerosol, CAD =100 reflects that the 
layer is most likely cloud, and CAD = 0 reflects inability to 
distinguish between aerosol or cloud (Liu et al., 2009). We 
required the CAD score to be less than —20 to reduce the 
possibility of including clouds misclassified as aerosol. 

In order to retrieve aerosol extinction from total attenuated 
backscatter, the quantity measured by the CALIPSO lidar, an 
initial lidar ratio (ratio of aerosol extinction to backscatter) 
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is assumed based on the aerosol type determined by the au- 
tomated aerosol subtyping algorithm (Young and Vaughan, 
2009; Omar et al., 2009). At times during the extinction re- 
trieval the initial lidar ratio may be adjusted to converge to 
a solution. An Extinction QC Flag is assigned for the ex- 
tinction retrieval of each aerosol layer which summarizes the 
final state of the retrieval solution. We required that Extinc- 
tion QC flag be 0, 1, 16, or 18, indicating that the either the 
initial lidar ratio was unchanged (0), the initial lidar ratio was 
measured based on the layer transmittance (1), the layer was 
opaque with an unchanged lidar ratio (16), or that the layer 
was opaque and the lidar ratio was reduced to avoid a di- 
vergent solution (18), respectively. Extinction uncertainty is 
reported along with each extinction value in the CALIPSO 
Level 2 aerosol profile product. A flag value of 99 000 Mm - 1 
is used to indicate that the extinction uncertainty solution is 
diverging to infinity, in which case the extinction retrieval 
should be discarded. We require extinction uncertainty to 
be less than 99 000 Mm -1 to remove these cases. Further- 
more, aerosol layers detected beneath clouds are not used be- 
cause cloud attenuation may significantly reduce the signal- 
to-noise ratio and impact the retrieved aerosol extinction so- 
lutions at lower altitudes. 

Acknowledgements. The authors wish to thank our pilot, Stan 
Unander, of the University of Illinois Institute for Aviation, for 
flying our profiles and being the on-site science technician, and 
Chuck Greenwood of The Greenwood Group for providing a safe 
and effective aircraft platform. We also thank Jim Wendell and 
Robert Albee of NOAA/ESRL for systems engineering and fabri- 
cation, and Derek Hageman of NOAA/ESRL for data acquisition 
and visualization software. We gratefully acknowledge funding for 
this project from the NOAA Climate Program Office. 

Edited by: A. Petzold 

References 

Anderson, T. L. and Ogren, J. A.: Determining aerosol radiative 
properties using the TSI 3563 integrating nephelometer. Aerosol 
Sci. Technol., 29, 57-69, 1998. 

Anderson, T. L., Covert, D. S., Wheeler, J. D., Harris, J. 
M., Perry. K. D., Trost, B. E., Jaffe, D. J., and Ogren, J. 
A.: Aerosol backscatter fraction and single-scattering albedo: 
Measured values and uncertainties at a coastal station in 
the Pacific Northwest, J. Geophys. Res., 104, 26793-26807, 
doi: 10. 1029/1999 JD900172, 1999. 

Anderson, T. L„ Masonis, S. J., Covert, D. S.. Charlson, R. J., and 
Rood, M. J.: In situ measurement of the aerosol extinction-to- 
backscatter ratio at a polluted continental site, J. Geophys. Res., 
105, 26907-26915, doi:10.1029/2000JD900400, 2000. 
Anderson, T. L., Charlson, R. J., Winker, D. M., Ogren, J. A., and 
Holmen, K.: Mesoscale Variations of Tropospheric Aerosols, J. 
Atmos. Sci., 60, 119-136, 2003. 

Andreae, M. O. and Gelencser, A.: Black carbon or brown car- 
bon? The nature of light-absorbing carbonaceous aerosols, At- 


mos. Chem. Phys., 6, 3131-3148, doi:10.5194/acp-6-3131-2006, 
2006. 

Andrews, E., Sheridan, P. J., Ogren, J. A., and Ferrare, R.: In situ 
aerosol profiles over the Southern Great Plains cloud and radia- 
tion testbed site: 1. Aerosol optical properties, J. Geophys. Res., 
109, D06208, doi: 10.1029/2003JD004025, 2004. 

Andrews, E., Sheridan, P. J., Fiebig, M., McComiskey, A., Ogren, J. 
A., Arnott, P., Covert, D., Elleman, R., Gasparini, R.. Collins, D., 
Jonsson, H., Schmid, B., and Wang, J.: Comparison of methods 
for deriving aerosol asymmetry parameter, J. Geophys. Res., Ill, 
D05S04, doi: 10.1 029/2004 JD005734, 2006. 

Andrews, E., Sheridan, P. J., and Ogren, J. A.: Seasonal dif- 
ferences in the vertical profiles of aerosol optical properties 
over rural Oklahoma, Atmos. Chem. Phys., 11, 10661-10676, 
doi: 1 0.5 1 94/acp- 1 1 - 1 066 1 -20 1 1 , 20 1 1 . 

Augustine, J. A., Hodges, G. B., Dutton, E. G., Michalsky, J. J., 
and Cornwall, C. R.: An aerosol optical depth climatology for 
NOAA's national surface radiation budget network (SURFRAD), 
J. Geophys. Res., 113, D11204, doi: 10.1029/2007JD009504, 
2008. 

Bergstrom, R. W.: Extinction and absorption coefficients of the at- 
mospheric aerosol as a function of particle size, Contr. Atmos. 
Phys., 46, 223-234, 1973. 

Bergstrom, R. W., Pilewskie, P., Pommier, J., Rabbette, M., Rus- 
sell, P. B., Schmid. B., Redemann, J„ Higurashi, A., Naka- 
jima. T., and Quinn, P. K.: Spectral absorption of solar radiation 
by aerosols during ACE-Asia, J. Geophys. Res., 109. D19S15, 
doi: 10. 1029/2003JD004467, 2004. 

Bergstrom, R. W., Pilewskie, P., Russell, P. B., Redemann. J., Bond, 
T. C., Quinn, P. K.. and Sierau, B.: Spectral absorption proper- 
ties of atmospheric aerosols, Atmos. Chem. Phys., 7, 5937-5943, 
doi: 10.5 194/acp-7-5937-2007, 2007. 

Berkowitz, C., Berg. L. K., Yu, X.-Y., Alexander, M. L, Laskin, A., 
Xie, Y., Jobson, B. T., Andrews. E.. and Ogren, J.: The influence 
of fog and airmass history on aerosol optical, physical and chem- 
ical properties at Pt. Reyes National Seashore, Atmos. Environ., 
45,2559-2568, 2011. 

Bodhaine, B. A.: Aerosol measurements at four back- 
ground sites, J. Geophys. Res., 88, 10753-10768, 

doi: 10. 1029/JC088iC15p 10753, 1983. 

Bodhaine, B. A.: Aerosol absorption measurements at Barrow, 
Mauna Loa and the South Pole, J. Geophys. Res., 100, 8967- 
8975, doi: 10. 1029/95JD005 13, 1995. 

Bodhaine, B. A. and Dutton, E. G.: A long term decrease in Arc- 
tic haze at Barrow, Alaska, Geophys. Res. Lett., 20, 947-950, 
doi:10.1029/93GL01 146, 1993. 

Bohren, C. F. and Huffman, D. R.: Absorption and Scattering of 
Light by Small Particles, John Wiley and Sons, New York, 530 
pp„ 1983. 

Bond, T. C. and Bergstrom, R. W.: Light absorption by carbona- 
ceous particles: An investigative review. Aerosol Sci. Technol., 
40, 27-67, doi: 10.1 080/0278682050042 1521, 2006. 

Bond, T. C., Anderson, T. L., and Campbell, D.: Calibration and 
intercomparison of filter-based measurements of visible light 
absorption by aerosols. Aerosol Sci. Technol., 30, 582-600, 
doi: 10. 1080/027868299304435, 1999. 

Charlson, R. J., Langner, J., Rodhe, H., Leovy, C. B., and Warren, 
S. G.: Perturbation of the northern hemisphere radiative balance 
by backscattering from anthropogenic sulfate aerosols, Tellus, 


Atmos. Chem. Phys., 12, 11695-11721, 2012 


www.atmos-chem-phys.net/12/11695/2012/ 



P. J. Sheridan et al.: Vertical profiles of aerosol optical properties over central Illinois 


11719 


43AB, 152-163, 1991. 

Clarke, A. D.. Shinozuka, Y., Kapustin, V. N., Howell, S., Hue- 
belt, B., Doherty, S., Anderson, T., Covert, D., Anderson, J., Hua, 
X., Moore II, K. G., McNaughton, C., Carmichael, G., and We- 
ber, R.: Size distributions and mixtures of dust and black carbon 
aerosol in Asian outflow: Physiochemistry and optical properties, 
J. Geophys. Res., 109, D15S09, doi:10.1029/2003JD004378, 
2004. 

Collaud Coen, M., Weingartner, E., Schaub, D„ Hueglin. C., Cor- 
rigan, C., Henning, S., Schwikowski, M., and Baltensperger, U.: 
Saharan dust events at the Jungfraujoch: detection by wavelength 
dependence of the single scattering albedo and first climatology 
analysis, Atmos. Chem. Phys., 4, 2465-2480, doi: 10.51 94/acp- 
4-2465-2004, 2004. 

Delene, D. J. and Ogren, J. A.: Variability of aerosol optical prop- 
erties at four North American surface monitoring sites, J. Atmos. 
Sci„ 59, 1135-1150, 2002. 

Dubovik, O. and King, M. D.: A flexible inversion algorithm 
for retrieval of aerosol optical properties from sun and sky 
radiance measurements, J. Geophys. Res., 105, 20673-20696, 
doi: 10. 1029/2000JD900282, 2000. 

Eck, T. F„ Holben, B. N„ Reid, J. S„ O'Neill, N. T„ Schafer, J. 
S., Dubovik, O., Smirnov, A., Yamasoe, M. A., and Artaxo, P: 
High aerosol optical depth biomass burning events: A compari- 
son of optical properties for different source regions, Geophys. 
Res. Lett., 30, 2035, doi:10.1029/2003GL017861, 2003. 

Eck, T. F., Holben, B. N., Dubovik, O., Smirnov, A., Goloub, P., 
Chen, H. B., Chatenet, B., Gomes, L., Zhang, X.-Y., Tsay, S.- 
C., Ji, Q., Giles, D. M„ and Slutsker, I.: Columnar aerosol op- 
tical properties at AERONET sites in central eastern Asia and 
aerosol transport to the tropical mid-Pacific, J. Geophys. Res., 
110, D06202, doi: 10. 1029/2004 JD005274, 2005. 

Ensor, D. S., Charlson, R. J., Ahlquist, N. C., Whitby, K. T., Husar, 
R. B., and Liu, B. Y. H.: Multiwavelength nephelometer mea- 
surements in Los Angeles smog aerosol - I. Comparison of cal- 
culated and measured light scattering, J. Colloid Interface Sci., 
39, 242-251, 1972. 

Esteve, A. R.. Ogren, J. A., Sheridan, P. J., Andrews, E., Holben, B. 
N., and Utrillas, M. P: Sources of discrepancy between aerosol 
optical depth obtained from AERONET and in-situ aircraft pro- 
files, Atmos. Chem. Phys., 12, 2987-3003, doi: 10.5 194/acp-12- 
2987-2012, 2012. 

Ferrare, R. A., Feingold, G., Ghan, S., Ogren, J., Schmid, B., 
Schwartz, S. E., and Sheridan, P: Preface to the special 
section: Atmospheric Radiation Measurement Program May 
2003 Intensive Operations Period examining aerosol proper- 
ties and radiative influences, J. Geophys. Res., Ill, D05S01, 
doi:10. 1029/2005JD006908, 2006. 

Fialho, P, Hansen, A. D. A., and Honrath, R. E.: Absorption coeffi- 
cients by aerosols in remote areas: A new approach to decouple 
dust and black carbon absorption coefficients using seven wave- 
length Aethalometer data, J. Aerosol Sci., 36, 267-282, 2005. 

Gasso, S., Hegg, D. A., Covert, D. S., Collins, D., Noone, K. 
J., Ostrom, E., Schmid, B., Russell, P. B., Livingston, J. M., 
Durkee. P. A., and Jonsson, H.: Influence of humidity on the 
aerosol scattering coefficient and its effect on the upwelling radi- 
ance during ACE-2, Tellus B, 52, 546-567, doi:10.1034/j,1600- 
0889.2000.00055.x, 2000. 


Hains, J. C., Taubman, B. F., Thompson, A. M., Stehr, J. W., 
Marufu. L. T., Doddridge, B. G., and Dickerson, R. R.: Origins 
of chemical pollution derived from Mid- Atlantic aircraft profiles 
using a clustering technique, Atmos. Environ., 42, 1727-1741, 
2008. 

Haywood, J. M. and Shine, K. P.: The effect of anthropogenic sul- 
fate and soot aerosol on the clear sky planetary radiation bud- 
get, Geophys. Res. Lett., 22, 603-606, doi:10.1029/95GL00075, 
1995. 

Holben, B. N., Eck, T. F., Slutsker, I., Tanre, D„ Buis, J. P., Set- 
zer. A., Vermote, E., Reagan, J. A., Kaufman, Y. J., Nakajima, 
T., Lavenu, F, Jankowiak, I., and Smirnov, A.: AERONET - A 
federated instrument network and data archive for aerosol char- 
acterization. Remote Sens. Environ., 66, 1-16, 1998. 

Holben, B. N., Tanre, D., Smirnov, A., Eck, T. F., Slutsker, 
I., Abuhassan, N., Newcomb, W., and Schafer, J. S.: An 
emerging ground-based aerosol climatology: Aerosol optical 
depth from AERONET, J. Geophys. Res., 106, 12067-12097, 
doi: 10.1 029/200 1 JD9000 1 4, 2001. 

Huebert, B. J., Howell, S. G., Covert, D., Bertram, T., Clarke, 
A., Anderson, J. R., Lafleur, B. G., Seebaugh, W. R., Wil- 
son, J. C., Gesler, D., Blomquist, B., and Fox, J.: PELTI: 
Measuring the passing efficiency of an airborne low tur- 
bulence aerosol inlet, Aerosol Sci. Technol., 38, 803-826, 
doi: 10. 1080/027868290500823, 2004. 

Intergovernmental Panel on Climate Change (IPCC): Climate 
Change 2007: The Scientific Basis, Contribution of Working 
Group I to the Fourth Assessment Report of the Intergovernmen- 
tal Panel on Climate Change, edited by: Solomon, S., Qin, D., 
Manning, M„ Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., 
and Miller. H. L., Cambridge Univ. Press, New York, NY, 2007. 

Kasten, F.: Visibility in the phase of pre-condensation, Tellus, 21, 
631-635, 1969. 

Kirchstetter, T. W., Novakov, T., and Hobbs, P: Evidence that 
the spectral dependence of light absoiption by aerosols is af- 
fected by organic carbon, J. Geophys. Res., 109, D21208, 
doi: 10. 1029/2004 JD004999, 2004. 

Koloutsou-Vakakis, S., Carrico, C. M., Kus, P., Rood, M. J., 
Li, Z., Shrestha, R., Ogren, J. A., Chow, J. C., and Wat- 
son, J. G.: Aerosol properties at a midlatitude Northern Hemi- 
sphere continental site, J. Geophys. Res., 106, 3019-3032, 
doi: 10.1 029/2000 JD900 1 26, 2001. 

Liu, Z., Vaughan, M. A., Winker, D. M., Kittaka, C., Kuehn, 
R. E., Getzewich, B. J., Trepte, C. R., and Hostetler, 
C. A.: The CALIPSO Lidar Cloud and Aerosol Dis- 
crimination: Version 2 Algorithm and Initial Assessment 
of Performance, J. Atmos. Ocean. Tech., 26, 1198-1213, 
doi: 10.1 175/2009 JTECHA1229.1, 2009. 

Marcq, S., Laj, P, Roger, J. C., Villani, P., Sellegri, K., Bona- 
soni, P, Marinoni, A., Cristofanelli, P., Verza, G. P., and Bergin, 
M.: Aerosol optical properties and radiative forcing in the 
high Himalaya based on measurements at the Nepal Climate 
Observatory-Pyramid site (5079ma.s.l.), Atmos. Chem. Phys., 
10, 5859-5872, doi: 10.5 194/acp-10-5859-2010, 2010. 

Massoli, P., Murphy, D. M., Lack, D. A., Baynard, T., Brock, C. A., 
and Lovejoy, E. R.: Uncertainty in light scattering measurements 
by TSI nephelometer: Results from laboratory studies and im- 
plications for ambient measurements. Aerosol Sci. Technol., 43, 
1064-1074, doi: 10. 1080/02786820903 156542, 2009. 


www.atmos-chem-phys.net/12/11695/2012/ 


Atmos. Chem. Phys., 12, 11695-11721, 2012 



11720 


P. J. Sheridan et al.: Vertical profiles of aerosol optical properties over central Illinois 


McFiggans, G., Artaxo, P., Baltensperger, U., Coe, H., Facchini. M. 

C. , Feingold, G., Fuzzi, S., Gysel. M., Laaksonen, A., Lohmann, 
U., Mentel, T. F., Murphy, D. M., O’Dowd, C. D., Snider, J. R., 
and Weingartner, E.: The effect of physical and chemical aerosol 
properties on warm cloud droplet activation, Atmos. Chem. 
Phys., 6 , 2593-2649, doi:10.5194/acp-6-2593-2006, 2006. 

McNaughton, C. S., Clarke, A. D., Howell, S. G., Pinkerton, M., 
Anderson, B., Thornhill, L., Hudgins, C., Winstead, E., Dibb, 
J. E., Scheuer, E., and Mating, H.: Results from the DC -8 In- 
let Characterization Experiment (DICE): Airborne versus surface 
sampling of mineral dust and sea salt aerosols, Aerosol Sci. Tech- 
nol., 41, 136-159, doi: 10. 1080/02786820601 118406, 2007. 

Moosmiiller, H., Chakrabarty, R. K., and Arnott, W. P: Aerosol light 
absorption and its measurement: A review, J. Quant. Spectrosc. 
Ra„ 110, 844-878, 2009. 

Muller, D., Wandinger, U., Althausen, D., and Fiebig, M.: Compre- 
hensive particle characterization from three-wavelength Raman- 
lidar observations: Case study, Appl. Optics, 40, 4863-4869, 
2001 . 

Muller, T., Nowak, A., Wiedensohler, A., Sheridan, P., Laborde, M., 
Covert, D. S., Marinoni, A., Imre, K., Henzing, B., Roger, J.-C., 
dos Santos, S. M., Wilhelm, R., Wang, Y.-Q., and de Leeuw, G.: 
Angular illumination and truncation of three different integrat- 
ing nephelometers: Implications for empirical, size-based correc- 
tions, Aerosol Sci. Technol., 43, 581-586, 2009. 

Myhre, G.: Consistency between satellite-derived and modeled es- 
timates of the direct aerosol effect, Science, 325, 187-190, 
doi: 10. 1 1 26/science. 1 1 7446 1 , 2009. 

Ogren, J. A.: In situ observations of aerosol properties, in: Aerosol 
Forcing of Climate, edited by: Charlson, R. J. and Heintzenberg, 
J., 215-226, J. Wiley. Hoboken, NJ, 1995. 

Ogren, J. A.: Comment on “Calibration and Intercompari- 
son of Filter-Based Measurements of Visible Light Ab- 
sorption by Aerosols”, Aerosol Sci. Technol., 44, 589-591, 
doi: 10. 1080/02786826.2010.4821 11, 2010. 

Omar, A., Winker, D., Kittaka, C., Vaughan, M„ Liu, Z., Hu, Y., 
Trepte, C., Rogers, R., Ferrate, R., Kuehn, R., and Hostetler, C.: 
The CALIPSO Automated Aerosol Classification and Lidar Ra- 
tio Selection Algorithm, J. Atmos. Ocean. Tech., 26, 1994-2014, 
doi: 1 0. 1 1 75/2009 JTECH A 1 23 1 . 1 , 2009 . 

Quinn, P. K., Bates, T. S., Baynard, T., Clarke, A. D., Onasch, T. 
B., Wang, W„ Rood, M. J., Andrews, E., Allan. J., Carrico, C. 
M., Coffman, D., and Worsnop, D.: Impact of particulate organic 
matter on the relative humidity dependence of light scattering: 
A simplified parameterization, Geophys. Res. Lett., 32, L22809, 
doi: 10. 1029/2005GL024322, 2005. 

Rahn, K. A., Borys, R. D., and Shaw, G. E.: The Asian 
source of Arctic haze bands, Nature, 268, 713-715, 
doi: 10. 1038/268713a0, 1977. 

Ramanathan, V., Crutzen, P. J., Lelieveld, J., Mitra, A. P, Althausen, 

D. , Anderson, J., Andreae, M. O., Cantrell, W., Cass, G. R., 
Chung, C. E„ Clarke, A. D„ Coakley, J. A., Collins, W. D., 
Conant, W. C., Dulac, F., Heintzenberg, J., Heymsfield, A. J., 
Holben, B., Howell, S., Hudson, J., Jayaraman, A., Kiehl, J. T., 
Krishnamurti, T. N„ Lubin, D, McFarquhar, G., Novakov, T., 
Ogren, J. A., Podgorny, I. A., Prather, K., Priestley, K., Prospero, 
J. M., Quinn, P K., Rajeev, K., Rasch, P., Rupert, S., Sadourny, 
R., Satheesh, S. K., Shaw, G. E.. Sheridan, P, and Valero, F. P 
J.: The Indian Ocean Experiment: An integrated assessment of 


the climate forcing and effects of the great Indo-Asian haze, J. 
Geophys. Res., 106, 28371-28399, doi: 10.1029/2001JD900133, 
2001 . 

Rogers, R. R., Hair, J. W., Hostetler, C. A., Ferrate, R. A., Obland, 
M. D., Cook, A. L., Harper, D. B., Burton, S. P, Shinozuka, Y., 
McNaughton, C. S., Clarke, A. D., Redemann, J., Russell, P. B., 
Livingston, J. M., and Kleinman, L. I.: NASA LaRC airborne 
high spectral resolution lidar aerosol measurements during MI- 
LAGRO: observations and validation, Atmos. Chem. Phys., 9, 
4811-4826, doi:10.5194/acp-9-481 1-2009, 2009. 

Rozwadowska, A.: Influence of aerosol vertical profile variability 
on retrievals of aerosol optical thickness front NOAA AVHRR 
measurements in the Baltic region, Oceanologia, 49, 165-184, 
2007. 

Schnaiter, M.. Gimmler, M., Llamas, I., Linke, C., Jager, C., and 
Mutschke, H.: Strong spectral dependence of light absorption 
by organic carbon particles formed by propane combustion, At- 
mos. Chem. Phys., 6, 2981-2990, doi:10.5194/acp-6-2981-2006, 
2006. 

Sheridan, P. J. and Ogren, J. A.: Observations of the vertical and re- 
gional variability of aerosol optical properties over central and 
eastern North America, J. Geophys. Res., 104, 16793-16805, 
doi: 1 0. 1 029/1 999 JD90024 1 , 1999. 

Sheridan, P. J., Delene, D. J., and Ogren, J. A.: Four years of con- 
tinuous surface aerosol measurements from the Department of 
Energy’s Atmospheric Radiation Measurement Program South- 
ern Great Plains Cloud and Radiation Testbed site, J. Geophys. 
Res., 106, 20735-20747, doi: 10. 1029/2001 JD000785, 2001. 

Sheridan, P. J., Jefferson, A., and Ogren. J. A.: Spatial variabil- 
ity of submicrometer aerosol radiative properties over the In- 
dian Ocean during INDOEX, J. Geophys. Res., 107, 8011, 
doi: 10.1 029/2000 JD000 166, 2002. 

Shinozuka, Y., Clarke, A. D., Howell, S. G., Kapustin, V. N., 
McNaughton, C. S., Zhou, J., and Anderson, B. E.: Aircraft 
profiles of aerosol microphysics and optical properties over 
North America: Aerosol optical depth and its association with 
PM 2.5 and water uptake, J. Geophys. Res., 112, D12S20, 
doi: 10. 1029/2006JD0079 1 8, 2007. 

Taubman, B. F., Hains, J. C., Thompson, A. M., Marufu, L. T.. 
Doddridge, B. G., Stehr, J. W., Piety, C. A., and Dickerson, 
R. R.: Aircraft vertical profiles of trace gas and aerosol pol- 
lution over the mid-Atlantic United States: Statistics and me- 
teorological cluster analysis, J. Geophys. Res., Ill, D10S07, 
doi: 10.1 029/2005 JD006 1 96, 2006. 

Turner, D. D., Ferrare, R. A., and Brasseur, L. A.: Aver- 
age aerosol extinction and water vapor profiles over the 
southern Great Plains, Geophys. Res. Lett., 28, 4441^1444, 
doi:10.1029/2001GL013691, 2001 . 

VanCuren, R. A.: Asian aerosols in North America: Extract- 
ing the chemical composition and mass concentration of the 
Asian continental aerosol plume from long-term aerosol records 
in the western United States. J. Geophys. Res., 108, 4623, 
doi: 10. 1029/2003JD003459, 2003. 

Vernier, J.-P, Thomason, L. W., Pommereau, J.-P, Bourassa, A., 
Pelon, J., Gamier, A., Hauchecorne, A., Blanot, L., Trepte, 
C., Degenstein. D., and Vargas. F.: Major influence of trop- 
ical volcanic eruptions on the stratospheric aerosol layer 
during the last decade, Geophys. Res. Lett., 38, L12807, 
doi: 10. 1029/20 11GL047563, 2011. 


Atmos. Chem. Phys., 12, 11695-11721, 2012 


www.atmos-chem-phys.net/12/11695/2012/ 



P. J. Sheridan et al.: Vertical profiles of aerosol optical properties over central Illinois 


11721 


Verver, G., Raes, F., Vogelezang, D., and Johnson, D.: The second 
Aerosol Characterization Experiment (ACE-2): Meteorological 
and chemical context, Tellus 52B, 126-140, doi:10.1034/j.l600- 
0889.2000.00090.x, 2000. 

Veselovskii, I., Kolgotin, A., Muller, D., and Whiteman, D. N.: In- 
formation content of multiwavelength lidar data with respect to 
microphysical particle properties derived from eigenvalue analy- 
sis, Appl. Optics, 44, 5292-5303, 2005. 

Virkkula, A., Backman, J., Aalto, R P., Hulkkonen, M., Riuttanen, 

L. , Nieminen, T., dal Maso. M., Sogacheva, L., de Leeuw, G., 
and Kulmala, M.: Seasonal cycle, size dependencies, and source 
analyses of aerosol optical properties at the SMEAR II mea- 
surement station in Hyytiala, Finland, Atmos. Chem. Phys., 11, 
4445-4468, doi:10.5194/acp-ll-4445-2011, 2011. 

Weingartner, E., Nyeki, S., and Baltensperger, U.: Seasonal and di- 
urnal variation of aerosol size distributions (10 < D < 750 nm) at 
a high-alpine site (Jungfraujoch 3580 masl), J. Geophys. Res., 
104, 26809-26820, doi:10.1029/1999JD900170, 1999. 

Winker, D. M., Vaughan, M. A.. Omar, A. H., Hu, Y., Pow- 
ell, K. A., Liu, Z„ Hunt, W. H„ and Young, S. A.: 
Overview of the CALIPSO Mission and CALIOP Data Pro- 
cessing Algorithms. J. Atmos. Oceanic Tech., 26, 2310-2323, 
doi: 1 0. 1 1 75/2009 JTECH A 1 28 1 . 1 , 2009 . 

Winker, D. M., Tackett, J. L., Getzewich. B. J., Liu, Z., Vaughan, 

M. A., and Rogers, R. R.: The global 3-D distribution of tro- 
pospheric aerosols as characterized by CALIOP. Atmos. Chem. 
Phys. Discuss., 12, 24847-24893, doi:10.5194/acpd-12-24847- 
2012, 2012. 


Wiscombe, W. J. and Grams, G. W.: The backscattered fraction 
in two-stream approximations, J. Atmos. Sci., 33, 2440-2451, 
1976. 

Young, S. A. and Vaughan, M. A.: The retrieval of pro- 
files of particulate extinction front Cloud Aerosol Lidar In- 
frared Pathfinder Satellite Observations (CALIPSO) data: Al- 
gorithm description, J. Atmos. Ocean. Tech., 26, 1105-1119, 
doi: 1 0. 1 1 75/2008 JTECHA 1221.1, 2009. 

Yu, H., Kaufman, Y. J., Chin, M., Feingold. G., Renter, L. A.. An- 
derson, T. L., Balkanski, Y., Bellouin, N.. Boucher, O., Christo- 
pher, S., DeCola, P, Kahn, R., Koch, D„ Loeb, N., Reddy, 
M. S., Schulz, M., Takemura, T., and Zhou, M.: A review 
of measurement-based assessments of the aerosol direct ra- 
diative effect and forcing, Atmos. Chent. Phys., 6, 613-666, 
doi:10.5194/acp-6-613-2006, 2006. 

Yu, H„ Chin, M„ Winker, D. M„ Omar, A. H„ Liu. Z„ Kittaka, C., 
and Diehl. T. : Global view of aerosol vertical distributions front 
CALIPSO lidar measurements and GOCART simulations: Re- 
gional and seasonal variations, J. Geophys. Res., 115, D00H30, 
doi: 10. 1029/2009JD0 13364. 2010. 


www.atmos-chem-phys.net/12/11695/2012/ 


Atmos. Chem. Phys., 12, 11695-11721, 2012 



